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Abstract 

The Hamilton-Lagrange action principle for Relativistic Schro- 
dinger Theory (RST) is converted to a variational principle (with con- 
straints) for the stationary bound states. The groundstate energy is 
the minimally possible value of the corresponding energy functional 
and the relativistic energy eigenvalue equations do appear as the cor- 
responding variational equations. The matter part of these eigenvalue 
equations is a relativistic generalization of the well-known Ritz prin- 
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ciple in non-relativistic quantum mechanics which however disregards 
the dynamical character of the particle interactions. If the latter are 
included in the proposed principle of minimal energy for the bound 
states, one obtains a closed dynamical system for both matter and 
gauge fields. The new variational principle enables the development 
of variational techniques for solving approximately the energy eigen- 
value equations. As a demonstration, the positronium groundstate is 
treated in great detail. Here a simple exponential trial function is suffi- 
cient in order to reproduce the (exact) result of conventional quantum 
mechanics where the relativistic and spin effects are neglected. 
PACS Numbers: 03.65.Pm - Relativistic Wave Equations; 
03.65.Ge - Solutions of Wave Equations: Bound States; 
03.65.Sq - Semiclassical Theories and Applications; 03.75.B 
- Matter Waves 
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I Introduction and Survey of Results 



The existence of an action principle is generally believed to be a very 
attractive feature of any (quantum) field theory, especially concerning the 
modern gauge field theories. The reason is that such an action principle does 
not only admit the deduction of the dynamical equations of the theory but it 
provides also many additional advantages, e.g. deduction of the conservation 
laws due to the symmetries of the theory (Noether theorem), canonical and 
path integral quantization, coupling to other fields, etc. (see, e.g., ref.s [HIS])- 
In view of such a considerable relevance of the variational principles, it does 
not appear as a surprise that there is an extended literature concerning the 
mathematical structure of the variational principles themselves, e.g. ref.s 
[31 i]. 

One important aspect of these variational principles refers to the possibil- 
ity of using them for the construction of certain approximation techniques if 
exact solutions of the dynamical equations cannot easily be found which will 
be mostly the case. The present paper is also concerned with just this aspect 
of the variational principles, namely in the context of the energy eigenvalue 
problems emerging within the framework of Relativistic Schrodinger Theory 
(RST), a recently estabhshed theory of relativistic quantum matter [5],[6l[7]. 
More precisely, the intention of the paper aims at the construction of a prin- 
ciple of minimal energy for the bound RST field configurations so that the 
groundstate carries the minimally possible value of the total field energy Et; 
and furthermore the excited states represent the stationary points {SEi[ = 0) 
of this energy functional i?T- It should be rather evident that the existence 
of such a minimal-energy principle is of invaluable practical usefulness for 
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the calculation of the (relativistic) energy levels of the bound systems (i.e. 
atoms and molecules). The point here is that one is not necessarily forced 
to look for the exact solutions of the eigenvalue equations but instead one 
can resort to the invention of appropriate trial configurations of the RST 
fields with all those symmetries (and other properties) being owned by the 
unknown exact solution, too. The principle of minimal energy can then be 
used in order to find the trial configurations with minimal (or stationary) 
energy which mostly is equivalent to a purely algebraic problem, namely the 
determination of the ansatz parameters in the chosen trial functions. 

For the sake of a simple demonstration and comparison with exactly solv- 
able examples of the conventional theory, we select the positronium ground- 
state as a typical two-body problem. Restricting ourselves here to the non- 
relativistic situation, one can easily show that the choice of a simple trial 
function reproduces exactly the positronium groundstate energy as it is pre- 
dicted by the conventional quantum mechanics. 

These results are elaborated through the following sequence of interme- 
diate steps: 

1. RST Eigenvalue Problem 

As the point of departure for the construction of the desired variational 
principle, one reconsiders the emergence of the typical energy eigenvalue 
problem within the framework of RST. Sect. II presents a brief sketch of 
the general theory for two oppositely charged particles such as, e.g., hydro- 
gen atom or positronium. The important point here is that the subsequently 
defined energy functional Et can be based upon the generally accepted hy- 
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pothesis of field energy wliich is concentrated in any relativistic field configu- 
ration with energy- momentum density T^^,, see equation flll.53|) below. This 
object Tf^i, may be deduced from the corresponding RST Lagrangean £rst 
pL57p in the usual way, i.e. via the standard Noether theorem. Indeed, the 
existence of a Hamilton-Lagrange action principle fill. 56ap - fill. 56bl) for RST 
is essential for the subsequent construction of the desired minimal-energy 
principle for the bound states since this action principle provides the possi- 
bility of introducing the Poisson identities which then work as constraints 
for the variational procedure (see the discussion of this point in the preceding 
paper [6]). 

Next, the stationary bound states are introduced in Sect. Ill and lead 
immediately to the mass eigenvalue equations fllll.lSp and fllll.lTp for the 
time-independent Dirac spinors ipa{'r) of both particles (a = 1,2). The inter- 
esting point with these mass eigenvalue equations refers to the fact that they 
can not only be deduced from the original Hamilton-Lagrange action principle 
flll.56al) - flll.56b[) by means of the stationary ansatz flIII.la|) - flIII.lbl) for both 
Dirac spinors, cf. flIII.14|) and flIII.161) . but these mass eigenvalue equations 
can also be deduced from an appropriately constructed mass functional M^c? 
flIII.24p . It is true, this mass functional can be understood to represent the 
immediate relativistic generalization of the well-known Hartree-Ritz vari- 
ational principle flIII.33aP - fjIII.33bp for factorized two-particle wave func- 
tions $(^1,^2) flIII.36p . but nevertheless Mtc^ cannot be accepted to repre- 
sent our wanted energy functional -Et because this mass functional Mpc^ suf- 
fers from the same deficiencies as the conventional Ritz principle; namely, in 
the latter approach the interaction between both particles is simply taken as 
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the rigid Coulomb interaction (see the conventional Hamiltonian (1III.34P ) 
whereas it is well-known that the gauge field, as the mediator of the relativis- 
tic interactions, must of course be treated as a dynamical object obeying its 
own field equations (see the Maxwell equations flll.24p ). 

And furthermore, there is a second deficiency inherent in the Hartree-Ritz 
approach which refers to the conventional dogma that wave functions should 
always be unique and non-singular. However, the RST treatment of positro- 
nium in the preceding paper ^ has shown that non-unique wave functions 
of the type flIII.601) must be admitted. Indeed, their use yields then a more 
realistic prediction of the positronium groundstate energy as compared to the 
Hartree-Ritz-Schrodinger approach, see ref. [6]. The conclusion from this is 
that in RST one must both admit a more general type of wave function and 
treat the particle interaction as a proper dynamical object. Otherwise one 
cannot expect to achieve a well-working principle of minimal energy which 
takes adequate account of both the matter and gauge fields! 

2. Exotic Quantum States 

The treatment of the positronium groundstate in the preceding paper [6] 
demonstrates that the minimal value of the RST energy functional can- 
not be reached by admitting exclusively these non-singular and unique wave 
functions as they are usually required by the Ritz-Hartree-Schrodinger ap- 
proach in conventional quantum mechanics. For instance, the requirement 
of physical equivalence of both positronium constituents (i.e. positron and 
electron) entails that any of the two fermions has vanishing spin component 
along the z-direction (equation pil.6ip ) which is quite unusual for fermionic 
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particles. As a consequence of this unusual behavior of the fermions, several 
other taboos of conventional quantum theory become broken, too: 

(i) The wave functions are singular at the origin (r = 0) and along the 
whole z-axis from the very beginning, see equations ( 1111. 62al) - ( 1111. 62b]) 



below; but nevertheless these singularities do not spoil the normaliza- 
tion conditions in the relativistic sense flll.44p . 

(ii) The wave functions become non-unique, e.g. in the sense of equa- 
tion (1111. 60P : but the observable physical densities (of charge, current. 



energy-momentum etc.) generated by these ambiguous wave functions 
are still unique and physically well-behaved. 

(ill) The interaction potentials, generated by the exotic states, are singular 
at the origin but less singular as the standard Coulomb potential so that 
their field energy is kept finite and can thus enter the wanted energy 



functional without causing infinities, see equations (llll.76p - (llll.77p 
below. 

(iv) The magnetic moment carried by the bound matter fields amounts only 
to half a Bohr magneton fiB, see the asymptotic form of the magnetic 
potential in equation (I1V.48|) below. 



3. Positronium Groundstate 

As a concrete demonstration, all this theoretical structure is evoked in 
order to calculate the positronium groundstate energy Eq. First, in the 
absence of an exact solution to the corresponding RST eigenvalue prob- 
lem (consisting of the coupled set of relativistic mass eigenvalue equations 



(lIII.66aP - (lIII.66dp and Poisson equations (lIII.47aP - (lIIL47d|) ) one resorts to 
a self-suggesting variational technique based upon the constructed energy 
functional Et fllV.lOl) which, however, is applied in this paper only in its 
non-relativistic approximation e!^^ flIV.26p . This means that one has to 
guess a trial function as realistic as possible (see the simple exponential wave 
amplitude R{r) (]V.5|) ) and substitutes this into the energy functional e!^^ 
(]IV.26|) . This energy functional is additively composed of the kinetic ener- 
gies i?kin(«) of both particles (lIV.20aP - (]IV.20bl) plus their electrostatic inter- 
action energy E^^^ flIV.22aD : the magnetic interaction energy Ei^^ flIV.22bD IS 
first neglected and treated afterwards as a small perturbation of the electric 
effects. 

The energy functional e!^^ f lIV.26p contains also two constraints which 
have to be respected for the deduction of the mass eigenvalue and Poisson 
equations as the variational equations due to that functional (i.e. Se!^^ = 0). 
The first constraint refers to the wave function normalization (as shown, e.g., 
by equation flIV.47P ) and is automatically satisfied by our ansatz fIV.Qap . 
However, the second constraint refers to the electric Poisson identities, such 
as fllll.SSp . and requires a more subtle argument: if one wishes to have the 
Poisson constraints also automatically satisfied by the trial functions one first 
has to solve the corresponding Poisson equations; i.e. equation (]III.72p for 
the present situation. But if all constraints are thus satisfied automatically 
by our trial function R{r), one substitutes this into the energy functional e!^'' 
(IIV.26P and obtains an ordinary function E!^\r^) of the ansatz parameter r^, 
i.e. equation flV.7p which according to the principle of minimal energy adopts 
the groundstate energy Eq fIV.Sp as its minimal value. This just coincides 
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with the corresponding prediction (IV. ip of conventional quantum mechanics. 
However, observe here that this groundstate energy Eq flV.l|) owns the sta- 
tus of exactness within the framework of the conventional theory, whereas 
in RST it appears as an approximation (even if all the relativistic effects 
including magnetism are disregarded) since our trial function R{r) fIV.Sp is 
surely not the exact solution of the non-relativistic RST eigenvalue problem 
in the electrostatic approximation. 

Finally, the magnetic interaction energy (1IV.56P is estimated in the 
lowest-order of approximation, equation flV.16P below. It turns out that 
the RST prediction for the hyperfine splitting of the positronium ground- 
states ^5*0 and ^5*1 amounts to only 1, 47- 10~^[ey], whereas the experimental 
value is 8, 41 ■ 10~^[el^] [H]. Thus this lowest-order RST prediction shows that 
for the hyperfine splitting one needs a better trial function. 
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II Two-Fermion Systems in RST 



In order to introduce the relevant notation, a brief sketch of the general 
two-particle theory is presented first so that the characteristic dynamical 
structure becomes obvious: matter dynamics, Hamiltonian dynamics, gauge 
field dynamics, action principle, and the associated conservation laws (for a 
more detailed presentation of the RST dynamics, see the preceding papers 
[S]-[7]). It is true, the existence of an action principle is common to almost 
all of the successful field theories, but a pleasant feature of the present RST 
dynamics refers to the fact that its action principle can be converted to a 
principle of minimal energy for the bound systems. This will subsequently 
be exploited in order to compute approximately the positronium groundstate 
energy. 

A. Matter Dynamics 

The central equation of motion for matter is the Relativistic Schrodinger 
Equation (RSE) 

ihcv^-^ = n^-^ , (II. 1) 

or if matter is to be described by an intensity matrix X in place of a pure 
state \1', one applies the Relativistic von Neumann Equation (RNE) 

v^i = ^^{in,-n,i) . (II.2) 

In the present paper, we will exclusively deal with pure two-particle states \E' 
which in RST are always the direct (Whitney) sum of the one-particle states 
ipa {a = 1,2), i.e. 

<il{x) =^i{x)®Mx) ■ (n.3) 
10 



Here the one-particle states ipai^) are four- component Dirac spinor fields so 
that the two-particle wave function \E'(x) may be understood as a section 
of a complex vector bundle over space-time as the base space with typical 
fibre C^. 

Both particles are interacting with each other via the principle of minimal 
coupling, i.e. the gauge- covariant derivative in the RSE (jll.ip is defined by 
means of the gauge potential {bundle connection) in the usual way as 

V^^ = d,^ + A,^ , (II.4) 

or, resp., in component form 

V^^ = (D^V-i) © {D,4^2) . (11.5) 

Here the gauge-covariant derivatives of the one-particle states ipaix) are given 
by 

Df,ipi = df^ipi - iAlipi - iB^i)2 (II.6a) 
D^ij2 = d^ij2 - iAliJ2 - iB^iJi , (H.6b) 

provided the bundle connection A^ takes its values in the four- dimensional 
Lie algebra U{2) of the unitary group U(2) {structure group) and is decom- 
posed with respect to a suitable basis of generators {tq, x, x} as follows: 

2 

A^ = Y, + B,x - B,x ■ (11.7) 

a=l 

The (real-valued) electromagnetic potentials do mediate the electromag- 
netic interactions between both particles; and similarly the (complex- valued) 
exchange potentials do mediate the exchange interactions which thus are 
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treated in RST as real forces on the same footing as their electromagnetic 
counterparts. However the exchange forces (due to 5^) can be active exclu- 
sively among identical particles and must vanish {B^ = 0) for non-identical 
particles (see refs.s O [7]). Since we restrict ourselves in the present paper 
to a system of two oppositely charged particles with different or identical 
masses Mp and Me, resp., the exchange forces must therefore be zero and 
consequently the covariant derivatives (lII.6aP - flll.6bp simplify to 



D^^i = d^^i - tAl^, (II.8a) 
D^ij2 = d^^2 - iAl^2 . (II.8b) 
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B. Hamiltonian Dynamics 



The Hamiltonian Ti^, occurring in the RSE OII.ip or in the RNE (111.2p . 
takes its values in the general linear algebra G£{2, C) and is itself a dynamical 
object which is to be determined from its field equations, i.e. the integrability 
condition 

v^n, - v,n^ + [n^, n,] = ihcj^^, (11.9) 

and the conservation equation 



-ihc 



Mc 



(11.10) 



The integrabihty condition flll.Qp contains the curvature T^;,v of the bundle 
connection flll.Tp 



(11.11) 



and guarantees the validity of the bundle identities 



(II.12a) 
(II.12b) 



The conservation equation OII.lOp contains the mass operator Ai and the 
Spin(l,3) generators S^,^ 

s^, = ^[r^,r,] (11.13) 

which both are assumed to be covariantly constant 



Vf^Yxu = 



(II.14a) 
(II.14b) 
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The latter constancy condition (111.14bp may be traced back to the covariant 
constancy of the total velocity operator 

VxT^ = , (11.15) 

where are the direct sum of the one-particle Dirac matrices 7^ 

r^ = (-7^)©7M (11.16) 

and therefore can be taken as the generators of the required eight-dimensional 
representation of the Clifford algebra C(l, 3), i.e. 

r^r, + r,r^ = 2g^, ■ i(8) . (ii.i?) 

Observe here that the arrangement of the plus and minus signs in the direct 
sum f lII.16p displays the opposition of both particle charges (positive charge 
of the first particle and negative charge of the second particle, by convention). 

The conservation equation (jll.lOp is needed for the deduction of the con- 
servation laws from the RST dynamics (see below) and admits an equivalent 
algebraic formulation: 

T^'Hf, = TY^r^ = Mc"^ . (11.18) 

This can be used in order to eliminate the Hamiltonian Ti^ by recasting the 
RSE (jll.ip into the two-particle Dirac Equation (DE) 

ihT^'V^,^ = Mc^ . (11.19) 

In component form, this equation reads 

ih-i^D^ipi = -MpCipi (II.20a) 
ih^^D^^2 = Me C^2 , (II.20b) 
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provided the (covariantly constant) mass operator A4 is written as 

2 

;w = i^MX (11.21) 

a=l 

with 

= Me (IL22a) 
= Mp , (IL22b) 

where Mp and Mg are denoting the rest mass of the positively and neg- 
atively charged particle, resp. For the case of pure states, one can elimi- 
nate the Hamiltonian also by differentiating once more the RSE (jll.ip 



and substituting therein the derivative of from the original conservation 
equation (lll.lOp which yields a second-order equation of the Klein- Gordon 
type (KGE): 

Vf'V^m + ( 1 ^ = -S^'^'-F^^^ . (11.23) 



However, subsequently we will prefer to deal with the first order equation flll.lQp 
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C. Gauge Field Dynamics 



In order to close the RST dynamics, one finally has to specify some 
field equation for the bundle connection Afj,. Our choice is the non-Abelian 
Maxwell equation 

V^J^^, = -Ama,J, (11.24) 

where the current operator may be thought to decompose with respect 
to the structure algebra basis {tq,, a = 1 . . .4} = {ra, x, x} follows: 

Jp^ = ir^^a = i V^i + i%^2 + gf,X - *9^X) ■ (11.25) 

Here, the Maxwell currents j"" ^{a = 1, 2) generate the electromagnetic poten- 
tials A"'^ (1IL7P which is seen by explicitly writing down the electromagnetic 
part of the general Maxwell equations 011.24^ in component form (a = 1,2) 



D'^FV=47r«sJ%. (11.26) 

The exchange currents j^^ = g^^ and j^^ = —g^ do generate the exchange 
potentials and (]II.7|) : but since we are dealing here exclusively with 



non-identical particles the exchange potentials B^^B^ must be put to zero 
so that the Maxwell equations (1II.26P become Abelian: 

V^F\, = A'Ka,j\ (IL27a) 
V^F\, = A'Ka,j\ . (IL27b) 

The formal reason for this is that the bundle curvature (]II.7p and its 
curvature J-'^y (III. lip 

2 

a=l 
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become projected onto the Abelian subalgebraZ//(l)©W(l) when the exchange 
fields B^, G^u are put to zero. 

D. Conservation Laws 

The right choice of the gauge field dynamics is not a trivial thing because 
it must be compatible with the already fixed matter dynamics (for both the 
pure states and the mixtures). However, this desired compatibility of our 
choice can be verified in the general case by the following arguments: First, 
the generally valid bundle identity 

VVJ^^^ = , (11.29) 

when applied to the Maxwell equations (111.240 . yields the following source 
equation for the current operator J'^ 

V^J^ = , (11.30) 

or in component form 

D^3\ ^ . (11.31) 

This means that the two-particle Maxwell currents j"" ^{a = 1, ... 4) must be 
constructed in terms of the two-particle wave function in such a way that 
the covariant source equations flll.31|) do actually hold just as a consequence 
of the RST dynamics! 

This compatibility requirement can be satisfied by first constructing the 
RST currents ja^ through 

3a^. = ^v^^m (11.32) 
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with the velocity operators v^^ being defined through the following anticom- 
mutators 

= 2 {^«' ^/^} • (11.33) 

Indeed, one can easily show that these RST currents ja/^ ( III. 321) do obey the 
source equations 

^^JaM = , (11.34) 

provided the wave function \1/ (or intensity matrix X, resp.) does satisfy 
the RSE ( III. II) (or the RNE (111.21) . resp.). However, observe here that the 
RST currents (III. 321) cannot a priori identified with the Maxwell cur- 
rents ^ (III.25P generating the gauge potentials via the Maxwell equa- 
tions (III.26p . Consequently, there must be established some link between the 
Maxwell currents ^ and RST currents ja^ in such a way that both source 
equations (]II.31|) and (111.340 are simultaneously valid! This requirement can 
be satisfied by conceiving j"" ^ and ja^ as contra- and covariant versions of 
one and the same object; namely by introducing a covariantly constant fibre 
metric Kap for the associated Lie algebra bundle 

DxK^p = , (11.35) 

and then putting 

Actually, such a compatibility tensor Kap can be found 

Ka(3 = Ci tr Ta " tr + C2 ti (r„ ■ Tp) , (11.37) 



(II.36a) 
(II.36b) 
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where Ci and C2 are constants which have to be chosen in such a way that 
the following constraint for the currents holds: 

2 2 

a=l a=l 

Here the total current of the two-particle system appears as the sum of the 
Maxwell (or RST) currents and acts as the source of the total electromagnetic 
field F^y 

F^u = + , (11.39) 

i.e. one easily deduces from the Abelian Maxwell equations (lll.27aP - flll.27bp 

the total Maxwell equation 

V^F^, = -47ra,j, . (11.40) 

Moreover, an immediate consequence of this Maxwell equation is the con- 
tinuity equation for the total current 

V% = . (11.41) 

Therefore the total charge z may be defined through 

j.dS^- , (11.42) 



'(5) 

which is independent of the chosen hypersurface {S) but must of course be 
zero because we are dealing with opposite charges. This may be realized more 
clearly by expressing the RST currents ja^ through the Dirac currents ka^ 

jifji = hfi = ^27m^2 (11.43a) 
j2/. = -kl^^ = -fpili^ipi ■ (11.43b) 
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Thus, since anyone of the two particles is assumed to carry just one charge 
unit, one will apply the following normalization of the wave functions for the 
stationary bound states (a = 1, 2) 

d^f (r) = 1 , (11.44) 

where the hypersurface [S) in fill. 421) is taken as a time slice {t = const.) of 
space-time; and the stationary form of the Dirac currents /c^^ (III.43aP - (III.43bp 
is of course 

ka,{x)= ((%(f);-4(r-)) . (11.45) 

Obviously the total charge z flll.42p becomes actually zero, namely by simply 
observing the sum requirement flll.38p and applying the normalization con- 
ditions (1II.44P together with the relationship fill. 43aP - (III. 43bp between the 



RST and Dirac currents. 

It is very instructive to consider the local charge conservation (111.4ip also 
from an other viewpoint: The total current may be defined alternatively 
through 

J,, = . (11.46) 

Carrying here out the differentiation process (111.4ip and using the Dirac equa- 
tion flll.lQp together with the covariant constancy of the total velocity oper- 



ator flII.15P actually yields just the total charge conservation law flll.4ip 



A similar procedure does apply also to the local energy-momentum conser- 
vation 

= , (11.47) 
where ^'^^T^,^ is the total energy-momentum density of the field configuration 
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and is composed of a matter part ^^''T^„ and a gauge field part ^'^^T^^ 

= (°)T,. + («)T,. . (11.48) 

The interesting point here is that the individual sources of the partial den- 
sities '^^^T^y and ^'^^T^^, turn out as the Lorentz forces which are mutually 
annihilating: 

(D)T^^ = _ = hcF^^yj^'' (11.49) 

so that the local law flll.47p can be true. However the crucial condition for 
this pleasant result is, that the RST dynamics (i.e. matter and gauge field 
dynamics) is chosen as described above and that the partial densities are 
defined as follows: 

(^)r,. = ^K^f^ [f%,F^' - ^g^^F\,F^->^^ (11.50a) 
(°)T^, = ^r^,^ , (II.50b) 

with the energy-momentum operator T^^, being given in terms of the Hamil- 
tonian and total velocity operator as 

T^, = -^{T^H, + n,T^ + T,H^ + n^T,) . (11.51) 

If matter can be described by a pure state \& (in place of a mixture) so that 
the matter density is given in terms of \E' by (1II.50b|) , then the Hamiltonian 
can again be eliminated from the matter density ^^''T^y by means of the 
DE fllLT9ll which yields 

ojl^f^ 

(11.52) 
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Clearly, the energy- momentum density ^^^T^^j is the crucial object for test- 
ing the practical usefulness of the theory, because the corresponding energy 
content i?T of the field configuration is given by the spatial integral of the 
time component *^^^Too(r), i.e. 

E^ = j d'f'^^^Tooif) . (11.53) 

Since the density fllL48|) appears as the sum of a matter and gauge 

field part, the same must hold also for the total energy E't (III.53P 

= + , (11.54) 

with the individual contributions being defined in a self-evident way as 

Ej, = j rfV(°)Too(r) (II.55a) 

Eg = j rf¥(G)Too(r-0 . (II.55b) 

Subsequently we will clarify the question whether for the groundstate of 
the stationary two-particle systems the energy functional flll.53p adopts its 
minimally possible value (-w principle of minimal energy). 

E. Action Principle 



The conservation laws for charge flll.4ip and energy-momentum fill. 471) 



can be directly deduced from the general RST dynamics, but a more elegant 
method is provided by the Noether theorem [9J. For the latter method one 
needs an action principle 

6Wkst = (I1.56a) 
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from which both the matter dynamics (1II.19P and the gauge field dynam- 
ics flll.241) may be deduced by the usual variational methods. The corre- 
sponding RST Lagrangean £rst has been specified as a sum of the matter 
part and gauge field part Cq f9] 

/:rst[^,^m] = + ^g[A^] , (11.57) 

where the matter part is given by 

tflc — — — 

Cj^m = — [^Tf" (P^^) - (P^^) r'^^] - ^A^c^^ (11.58) 

and the gauge field part by 

he 

Cg[A,] = K^pF^^F^^-" . (11.59) 

Concerning the latter part ( 111.59^ . observe here that the bundle curvature 
(III. lip takes its values in the subalgebra W(l) ©W(l) because we are 
dealing with non-identical particles; and thus the gauge field Lagrangean 
becomes reduced to 

2 



fir 

Cg[A,] = V K,,F\,F'^^ . (11.60) 

a,b=l 



If the self-interactions are neglected, the fibre submetric Kab is of a very 
simple shape [5] 

{Kab} =(^^^ j . (11.61) 

Thus, from the formal point of view, /^Gl-^Ai] describes the interaction of the 
two gauge field modes F"^,^ (a, 6 = 1, 2), i.e. 

Cg[A^] = ^ [e, 'E2-H1' H2) , (11.62) 
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provided the four-tensors -F"^,^ are splitted into their space and time compo- 
nents Ea, Ha as usual 

Ea = {^''^E^} = {F\} (II.63a) 
Ha = {^''^W} = l^e^^F'^fc'l . (II.63b) 

But once the Lagrangean has been specified, it is an easy exercise to 
deduce both the RST matter dynamics flll.lQP and the gauge field dynam- 
ics f lll.24p from the action principle flll.56ap - flll.56bl) as the corresponding 
Euler-Lagrange equations. Furthermore, the considered conservation laws of 
charge fill. 300 and energy-momentum f]11.47p are just those which are pre- 
dicted by the Noether formalism, see ref. [9]. 
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Ill Stationary Bound Systems 



In order to make the proposed exercise with the Euler-Lagrange equations 
somewhat more instructive and fruitful, one may immediately pass over to 
the stationary systems which are defined through the usual product ansatz 
for the wave functions 

il>i{f,t) = exp [—i — — 1 \ ■ il)i{r) (Ill.la) 

^2(r,t) = exp (-^^^) ■ V'2(r) , (IILlb) 

whereas the electromagnetic potentials become time- independent: 

A\=[^'^^A,{r)--Aa{r)] (III.2) 
(a = 1,2) . 

This time-independence obviously does then apply also for the Dirac cur- 
rents kaf, flII.43ap - flII.43bD . see equation fllL45|l . 

A. Mass Eigenvalue Equations 

The mass eigenvalues Ma {a = 1,2), occurring in the stationary ansatz 
pil.lall -f llll.lbj) . must be determined through solving the stationary form of 



the matter dynamics fill. 20a|) - fill. 20b|) . This stationary form {^^mass eigen- 



value equations") may be obtained either by direct substitution of the station- 
ary ansatz flIII.laP - f]III.lb|) into the coupled Dirac equations flll.20ap - fjll.20bp 



or by substitution of that ansatz into the matter Lagrangean £d[^] flll.SSP 
and then carrying out the variational procedure with respect to the spatial 
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parts ipai'f') of the wave functions. Here it is easy to see that the mat- 
ter Lagrangean £d[^] sphts up into a sum of the single-particle contribu- 
tions £0(0), i-e. 

Cj,m = CM + Cj,{2) , (111.3) 
with the individual contributions being given by 

CM = -thc^iYiD^^,) - Mp c^V^i^i (111.4a) 
CM = ihci,2l'\D^ij2) - c2V>2^2 . (ni.4b) 

However, observe here that these matter contributions to the Lagrangean 
-^RST do nevertheless contain the gauge fields which invade the matter La- 
grangean via the gauge- covariant derivatives D^ipa (lll.8aP - (111.8bp . There- 



fore the matter contributions themselves split up into the proper kinetic 
parts (a) and the electric plus magnetic interaction contributions in the 
following way: 

CM = 4'"V) + + 4"V) + 4""^) , (ni.5) 

i.e. for the first (positively charged) particle 

C'i^\^^ = -ihc^i(f)^'V^i{f) (in.6a) 

4\i) = ■ ^'^^o(r) (II1.6b) 

C^^\^^= hcA2{f)-ki{^ (111.6c) 

£Sf = -Mic^ ■ (i)fco(r-) - Mpc^ ■ V^i(rl^i(rl (111.6d) 
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and similarly for the second (negatively charged) particle 



-^D (2) 



he (^ko(f) ■ ^%^{r) 



(m) 



D (2) - -he Ai{r) ' k2{r) 

= M2C^ ■ (%(f) - Me ■ ^2(r)V'2(r) 



^(M) 
-^D (2) 



(II1.7a) 
(IIL7b) 
(III.7c) 
(lll.Td) 



(III.8) 



From this reason, the space part (Wd, say) 

of the matter contribution to the action integral Wrst ( 1II.56b[) is not only 
built up by the kinetic and rest mass terms but contains also the electric (e) 
and magnetic (m) mass equivalents Mjjj^-'c^ of the gauge field energy, i.e. for 
the first particle 

(111.9) 



and similarly for the second particle 

-(2) 



(111.10) 

Here the mass renormalization factors Z'^^^ are defined through (a = 1, 2) 

rf¥V',(f)7/',(f) . (111.11) 

Furthermore the electric mass equivalents of the interaction energy appear 



cr2 



as 



M'f^c^ = hcj d^f (%o(r) ■ '^^^{r) 
Mifc^ = -hc I d^r (%o(r) ■ ^^\{r) 



(111.12a) 
(111.12b) 
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and analogously for the magnetic mass equivalents 

^M^2 _ j £fX2{r) • fci(f) (III.13a) 
Mjf ^c^ = ficj c/¥li(f) • A?2(r) • (111.13b) 



But with this explicit structure of the matter Lagrangean (]III.3|) it 



is easy to see that the variation of the action integral VTrst (lII-56aP -( ]II.56bp 
with respect to the first wave function ipi 

5(l)W^RST = 5(1)W^D(1) = (III.14) 

yields just the first mass eigenvalue equation 

ij • Vtpi (f) + (% (r)7>i (f) - ^2 (r) • (r ) = - 7o + ■ (r) , 

(III.15) 

and similarly the variation of Wrst with respect to the second wave func- 
tion iIj2 

° I 

S(2)WrsT = ^(2)^^0(2) = (111.16) 

yields the second mass eigenvalue equation 

n'ViJ2{r) + (%(f)7V2(0 - ^i(^) • 7^2 (^) = - (^—^ 7o - -y-^ ■ ^2(r) . 

(III.17) 

B. Mass Functional 

An interesting property of the mass eigenvalue equations (1111.150 and 
( IIII.17I) refers to the fact that they are linear with respect to the wave func- 
tions tpa- This implies that one can multiply any solution tpa{r) by some 
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constant and then obtains a further solution, which necessarily can not mod- 

o 

ify the value of the matter functionals VFD(ct) upon the whole class of such 

o 

solutions. On the other hand, the matter functionals WY){a) flIII.8|) are bilin- 



ear with respect to the wave functions and their (pseudo-) Hermitian conju- 

_ o 

gates tpa] and from this one concludes that Wj){a) must take the value zero 
upon the solutions of the mass eigenvalue equations (a = 1, 2): 

WM = . (III.18) 

This is the reason why one is not forced to impose some normalization con- 
dition upon the wave functions ipai"^ when one deduces the mass eigenvalue 
equations as the Euler-Lagrange equations due to the RST variational prin- 
ciple flII.56aD - flII.56bD . 

o 

However, it is just this property flIII.18|) of the matter functionals W]:,{a) 
which enables one to resolve these equations for the mass eigenvalues Ma, 
where it is convenient to apply the normalization conditions f ill. 441) for the 
wave functions i'air); and these mass functionals ( M[a]C^, say) emerge then 
in the following form: 

-M[i]c2 = ) ■ Mp + 2Tkin(i) + Mj^'^c^ + m/'^^c^ (III.19a) 
M[2]C2 = Zf2) ■ Me + 2Tkin(2) + M^^^c^ + ^c^ . (111.19b) 

Here the mass renormalization factors Z"^^-^ have already been specified by 
equation fllll.lip and the mass equivalents of the interaction energies by 
f lIII.12aP - fjIII.13bp . The emergence of a pre-f actor of two in front of the 
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kinetic energies Tkin(o.), being given by 

71.in(i)= '-hcjd^fij^^-VM^ (III.20a) 

Tkin(2) = -^^c j d^f^27'^M^ , (III.20b) 

is a relativistic effect and is compensated by the mass renormalization fac- 
tors Z'^^^ (see the discussion of this effect in ref.s [Sllin])- But clearly, in place 
of resolving the equations (IIII.ISP for the mass eigenvalues Ma in order to 
obtain the mass functionals Mf^ic^ (IIIL19aP - flIII.19bp . one could equally well 



multiply through the original mass eigenvalue equations ( llll.lSj) and ( lIII.lTj) 
by ipi and '?/^2 and integrating over whole three-space in order to arrive again 
at the same mass functionals M[a]C^ (IIII.19aP - (IIII.19b|) . 



The relevance of the mass functionals originates now from the fact that 
they are stationary upon the solutions of the mass eigenvalue equations (lIII.lSp 
and (lIII.17p . albeit with regard of the constraints of wave function normal- 
ization (III.44p . In order to take account of these constraints, one introduces 



the Lagrangean multipliers Ad (a) and considers the modified mass function- 
als M[a]C^ 

-M[i]c2 = -M[i]c2 + Ad(i) ■ ArD(i) (III.21a) 

M[2]C' = M[2]C2 + Ad(2) ■ A^D(2) , (III.21b) 

with the constraints Nj){a) being given by the wave function normaliza- 
tion (a = 1, 2), i.e. 

NM ^ J d^riJal'^^IJa^^ -1 = 0. (III.22) 

Indeed, carrying through here the variational procedure and comparing the 
emerging variational equations to the former mass eigenvalue equations (IIII.lSp 
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and (1III.17P just yields the identity of both sets of equations, provided one 
fixes the Lagrangean multiphers Ad (a) in terms of the mass eigenvalues Ma 
as follows: 

AD(y = Mic^ (IIL23a) 
Ad(2) = -Mac^ . (III.23b) 

Thus the Lagrangean multipliers just turn out to be identical to the mass 
eigenvalues (up to sign). Clearly such a result strongly remembers one of 
the well-known Ritz variational principle of conventional quantum mechan- 
ics which has frequently been applied in the early days of atomic physics [11] 
and in the meantime has been advanced to a standard exercise for any stu- 
dent [12]. 

Indeed, anyone of the mass functionals M\a]C^ flIII.19a|) - flIII.19b|) appears 



to be of a very plausible form, namely as the sum of three kinds of energy: 
rest mass energy, kinetic energy and interaction energy of the electric (e) 
and magnetic (m) type. Therefore it is very tempting to think that the 
total energy of the bound two-particle system could be identified with the 
sum (Mtc^) of both mass eigenvalues 

M^c^ = -M[i]c2 + M[2]c2 , (III.24) 

especially because the corresponding variational equations just coincide with 
the mass eigenvalue equations (IIII.ISP and (lIII.lTp . However, such a supposi- 
tion is incorrect; and a closer inspection of how the Ritz variational principle 
emerges as the non-relativistic limit of the present mass functional approach 
will reveal the origin of the deficiencies connected with the non-relativistic 
Ritz method. Furthermore this analysis provides one with a valuable hint 
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on how to overcome the shortcomings of those non-relativistic variational 
methods and to construct the relativistic energy functional Ei. 

C. Ritz Variational Principle 

In order to deduce the non-relativistic limit form of both the mass func- 
tionals Mj^jc^ fllll.19ap - fllll.19bl) and of their variational equations fllll.lSp 
and f]111.17p one conceives the Dirac four-spinors ipai'f^ as a direct sum of 
two-component Pauli spinors *^"V±(^) {(^ = 1;2), i.e. one puts 



and then one deduces the corresponding eigenvalue equations for these Pauli 
spinors from the original mass eigenvalue equations fllll.lSp and fllll.l7p . This 
yields for the first particle (a = 1) [10] 



- (a) 



(111.25) 



Mp - M] 

h 



(111.26a) 

-c-(^V+(r) 



h 



(111.26b) 



and similarly for the second particle (a = 2) 



la-V^^'^ip+if) + (%(f) ■ (2)y,_(f) - Ai • a('V+(^) 



Af c + Ma 

h 



c 




(111.27a) 



Me - Ma 



c 



■ (^V-f (r) . 



(111.27b) 
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Next one approximately resolves the eigenvalue equations (lIII.26aP and 
(IIII.27aP for the "small" components '•''V-(r) 



«^_(r)^ ^a.V(i)^+(r) (IIL28a) 

^ -^^'^^'^Mrl (ni.28b) 

2Mp c 



and substitutes this into the eigenvalue equations ( ]III.26bl) and ( ]III.27b[) for 
the "small" components (r) in order to obtain the well-known Pauli 

equations for the "large" components: 

-^A(i)^+(f) + hc^'Kir-} ■ (^V+(r1 = i?s(i) ■ ^'V+(r1 (III.29a) 

-^A(2)y,4f) - ;ic(%(f) ■ (^)^+(f) = Es(2) ■ ^'V+(r) . (III.29b) 

Here the Pauli-Schrodinger energy eigenvalues Es{a) are introduced through 

Es(i) = - (Mp + Mic^) (111.30a) 
Es(2) = M2c2-Mec2 , (IIL30b) 

which again demonstrates that the mass eigenvalue of the first particle is neg- 
ative (Ml < 0), in contrast to the second eigenvalue (M2 > 0). Furthermore 
the magnetic interaction terms (~ Aa) are omitted since the corresponding 
magnetic interaction energy is mostly much smaller than its electric coun- 
terpart described by ^"•^Ao^r). As a consequence of this omission, the Pauli 
spinors can be assumed to occupy a fixed direction in spinor space, e.g. 

('^)^+(r) = ^,(r) ■ I T> , (111.31) 

so that for the scalar wave functions (pa{r) there emerge the conventional 
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Schrodinger equations from the Pauli system (1111.29ap - (1111.29bp : 

-^A^i(f) + hc^'koif) ■ ^i(f) = i?s(i) ■ ^i(rl (III.32a) 



A^r) - hc^'koif) ■ if^if) = Es(2) ■ <f2{r) . (III.32b) 



2Me 

Observe here that either of the two particles feels the attractive force due to 
the other one because the first potential (^ko(r) (due to the positively charged 
particle) is positive and the second potential *^^ko(r) is negative! (See below 
for the discussion of the corresponding Poisson equations (lIII.47aP - (lIII.47dp ). 

Now the interesting point with this non-relativistic approximation (]III.32a|) - 
(1III.32bl) to the properly relativistic eigenvalue equations ( 1III.15I) and ( 1III.17I) 
refers to the fact that the non-relativistic system may be also deduced from 
a variational principle; i.e. the well-known Ritz principle |T3] 



SWs = (III.33a) 
Ws = [[ d'n d^r2 l>(ri, f2)Hs^ri, fs) , (III.33b) 



where the conventional Hamiltonian Hs is given by 



Hs = ^ + ^- . (III.34) 

^ 2Mp 2Me ||ri - rail ^ ' 



Indeed, carrying out the variational procedure (]III.33ap . with the constraint 



of wave normalization, lets emerge the conventional Schrodinger equation as 
the corresponding variational equation 

H^Wi.r^) = Es ■ $(fi,f2) . (111.35) 

It is generally believed that the associated conventional eigenvalue Es is the 
"true" (albeit non-relativistic) energy of such a two-particle system being 
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specified by the Schrodinger Hamiltonian (1111.340 ; and consequently all other 
predictions can at most appear to be approximations to that true value! For 
instance, one may try (for the groundstate) the following product ansatz 
(Hartree approximation) 

<l>(fi,f2) = (/^i(fi)-^2(r2) (III.36) 

and use this for carrying out the variational procedure (lIII.33ap - (]III.33bp 
which then yields the following one-particle eigenvalue equations: 

-^A^i(f) + he (2)\/Hs(r-) ■ ^i(r-) = -As(i) ■ (III.37a) 
-A(^2(rl - he (^VHs(rO ■ <^2(r) = -As(2) ■ V2{^ ■ (III.37b) 



2Me 

Here the normalization conditions for the conventional wave functions V5a(^) 

j d^r^a{^^a{^-l = Q (111.38) 

are respected by application of the method of Lagrangean multipliers; and 
the Hartree- Schrodinger potentials ^"■'VHs(r) given in terms of the one- 
particle wave functions ^ai'f^) through 

")V„.M= aJSr^'^itllpl („,.39a) 

f2)T/ /-N /"j3- ^^2(^2) ■ V52(^2) 



/ks(r) = -«s / d'r2 Z: . ^' . (III.39b) 

Clearly, the Lagrangean multipliers As (a) in the Hartree-Schrodinger eigen- 
value equations (lIII.37ap - (]III.37b|) are to be identified again with the energy 
eigenvalues Es{a) (a = 1,2) 

As(a) = -Es{a) , (III.40) 
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cf. the relativistic version (lIII.23aP - (lIII.23bp : and thus the Hartree-Schrdding 



system flIII.37a|) - flIII.37b|) is identical to the non-relativistic RST limit UII. 32a\) - 



mi. 32b\) . provided one can show (see below) that the electric RST poten- 
tials ("Uo(r) are identical to the Hartree-Schrodinger potentials ^"'^Vusif^ 
f lIII.39ap - flIII.39bl) . This important identification of the Hartree-Schrodinger 
(or more generally: Hartree-Fock) approach with the non-relativistic limit of 
the RST eigenvalue equations suggests that the RST variational method 

5Mt = (III.41) 

(with the mass functional Mt given by equations f]III.21aP - f]III.24p ) may be 
considered as a viable relativistic generalization of the Hartree-Ritz varia- 
tional principle (or Hartree-Fock approach, resp). But since the latter ap- 
proach is in general used as merely an approximation to the conventional 
Schrodinger theory, being based upon the "exact" non-relativistic eigenvalue 
equation pH.35p . it may seem that RST is afflicted with all those deficien- 
cies of the Hartree-Fock approach, albeit on a relativistic level (for a critical 
discussion of the limits of the HF approach see ref. [H]). However this 
conclusion is not valid because RST is capable of overcoming the main de- 
ficiencies of both the Ritz principle and the HF approach (apart from their 
non-relativistic character). These deficiencies are the following: 

i) the Ritz principle f]III.33aP - f]III.34p relies upon the instantaneous Cou- 
lomb interaction of the particles and thus violates the true spirit of 
relativity 

ii) the Hartree-Fock approach relies upon unique wave functions ^Pai^^ 
pn.36p . whereas not the wave functions themselves (as unobservable 
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objects) need be unique but rather the observable quantities (i.e. the 
physical densities in RST). 

In the following, we will explicitly demonstrate that the treatment of the 
interaction fields as proper dynamical variables together with the use of non- 
unique wave functions may actually lead to predictions which can compete 
with those of the conventional quantum theory. Indeed, this improvement 
of the standard Ritz-Hartree-Fock approaches will then result in the desired 
principle of minimal energy for the stationary bound states. 
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D. Poisson Equations 



Surely, if the particle interaction is to be considered as a dynamical object, 
then it must obey some equation of motion and cannot be specified by the 
rigid Coulomb interaction. However in RST, the existence of such a (relativis- 
tic) equation of motion for the interaction fields is a matter of course, since 
the RST variational principle (]II.56aP - (lII.56bp includes the gauge field Afj, a 
priori. It is merely necessary to carry out the variational procedure (111.56ap 
with respect to the bundle connection which then yields the general 
Maxwell equations (111.240 . or their Abelian specialization (lll.27ap - flll.27bp . 
resp. Since for the present Abelian situation (for non-identical particles) the 
field strengths F""^^ degenerate to the simple curls of the potential ^ 

F\, = d^A% - d,A\ (111.42) 

(see ref. |^ for the case of identical particles) the Abelian Maxwell equa- 
tions flll.27ap - flll.27bp yield just the well-known wave equations of classical 
electrodynamics (a = 1,2): 

d^d^A\ = A'Ka,j% , (III.43) 

provided one applies the usual Lorentz gauge condition 

(9^A% = . (III.44) 

For the present stationary situation flIII.2p . this splits up into the (electric) 
time component (a = 1,2) 

A(''Uo(r-) = -47ra, ^%{r) (III.45) 



38 



and (magnetic) space component 

AX(r) = -AnaJai^ ■ (ni.46) 

Or, if the Maxwell currents = {^'^^jo{r), —jair)} are written in terms of 
the stationary Dirac currents fca^ (111.45^ one has the Poisson equations 

A(^ko(r) = -Ana, '■%{r} (III.47a) 

A(2ko(f) = 47ras ^^^fco(f) (III.47b) 

AAi(f) = -47r«s A?i(f) (III.47c) 

AA2{r} = Aira, k2{f) . (III.47d) 



Observe here that the fibre submetric Kab 0II.37P has just the right form flll.6ip 



in order to get the Poisson equations adapted to the positive and negative 
charges carried by the particles! The standard solutions of the Poisson equa- 
tions are adopted to be of the usual form 



(iko(f) = a, / d^f' ^'^Vn (III.48a) 
(%(f) = -a, [ (IIL48b) 



A^(r) = a, I d^f'J^^i^ (IIL48c) 
/2(f) = -as j rfV^^^ . (IIL48d) 

The non-relativistic limits of these formally exact solutions are obtained by 
substituting herein the non-relativistic approximations for the Dirac densities 
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of charge ^""^ko (r) and current ka{r) [7]: 

(f) = MrhoMr) = ^"VU^) • ^"V+(^) + • ^"V-(^) (in.49a) 

(III.49b) 

Observe here that for the non-relativistic hmit the "small" Pauli com- 
ponents ^"^(y9_(r) can be neglected against its "large" counterparts ("^y9+(r*) 
only for the charge density ^""^ko^r) flIII.49ap . but not for the current den- 
sities flIIL49bp as the sources of the magnetic fields! This is the reason 



why we omitted the magnetic terms when deducing the Hartree-Schrodinger 
eigenvalue equations (]III.32aP - (]III.32bp for the "large" components from their 
properly relativistic RST form (jlll.lSp and (1111.170 . But inserting now the 
approximate form fllll.49ap of the charge densities ^""^o (r) into the formal 



solutions ^"'^Ao{r) flIII.48aP - flIII.48bp of the Poisson equations yields 



(%(f) ^ -a. I d^f' ^^^p^Jp^p , (111.50b) 

and this is just the result which is needed in order to identify the non- 
relativistic RST limit fjIII.32aP - fjIII.32bp with the conventional Hartree-Schrodinger 
eigenvalue system flIII.37aP -f lIII.37bp . 

This is the way in which RST cures the first one (i) of the Ritz-Hartree- 
Fock deficiencies mentioned above, namely by adopting the Poisson equations 
(or more generally: the Maxwell equations) for the determination of the 
gauge potentials and thus accepting the action of the gauge forces as a truly 
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dynamical process. However, observe in this context that the Poisson equa- 
tions emerge here as the variational equations due to the Hamilton-Lagrange 
action principle flll.56al) - flll.56bl) . not due to the variation of the RST mass 
functional Mt(? f lIII.24p ! Therefore it is not yet possible to identify this 
mass functional with the wanted energy functional as the ultimate goal of 
the present investigation. This goal will be readily achieved in the next sec- 
tion but can be prepared here by first regarding an important property of 

o 

the time-independent gauge part Wcl-^Ai] 

Wg[A^] = j d?rCG[A^] (HI.51) 

of the original RST action principle (]II.56bp . Indeed, substituting here the 
Lagrangean density Cq in the form flll.62p with the electromagnetic three- 
vector fields Ea{r) and Ha{r) being written in terms of the corresponding 
potentials ^""^Ao^r) , Aa{r) as 

Ea{r) = -V("Uo(f) (HI. 52a) 

Ha{f) = V X Aa{f) (III.52b) 

o 

lets appear that time- independent gauge part Wcf-^/i] the following form: 
Wc[A,] = rfV{v(%(f). V(%(f) - (V X A(f)) . (V X l2(r))} 

(ni.53) 

Furthermore, by resorting to the gauge field contributions C^^'^^ (lIII.6bp . 
(]III.6cP and (]III.7bp . (lIII.7cP due to the matter Lagrangean £d, one has the 
corresponding electric (e) and magnetic (m) action constituents as 

W^^ = -he [ d^f {(%(f) ■ (i)A;o(r) - (%(f) ■ (%(f)} (III.54a) 



= he d^r lA2{r)'ki{r) - Ai{f)'k2{r)\ . (lllMh) 
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Thus the former Poisson equations (lIII.47ap - (lIII.47d|) are actually recovered 

o o (e) o (m) 

from here by variation of the partial sum of action integrals Wq + Wj^ + Wj^ 
with respect to the static gauge potentials. 

Now the interesting point with this variational procedure for the gauge 
fields is that it leads us to global identities which subsequently will be needed 
as constraints for the principle of minimal energy, i.e. the Poisson identities. 
These global relations between the gauge fields and their sources emerge from 
the Hamiltonian-Lagrange action principle (lII.56aP - (]II.56bp by considering 
the scaling variations for the potentials, e.g. for the first electrostatic potential 

(%(r-) -> ^%{r) = a ■ (%(f) , (III.55) 

with the scaling factor being a constant over three-space. Similar ar- 
guments do hold also for the other potentials '^^]4o(r), v4a(r). By inserting 
this special variation (1111.55^ into the RST action integral, its relevant parts 
become 



o (e) o (e) 



he 



Ana 



j rf¥ {v (^ko(r) ■ V (%(f) + At: a, ^^%{r) ■ ^%^{ 



. (III.56) 



Since the Hamiltonian-Lagrange action principle demands stationarity of the 
action integral VFrst with respect to the choice of C^,: 



da 



c*=i 



, (III.57) 



one concludes from equation ( jIII.561) that the following integral relation must 
hold: 



iV(f)(i) = / d'f 



V(%(f)) • fv(%(f)) +47ras^%(r)-(%(f 



= , 
III.58) 
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and analogously for the other gauge potentials 
N^^\2) = [ £r [(v(%(f)) • (v(%(f)) -47ras(%(f) ■(%(f) 



n''^\2) = I d^f ( V X Ai(f) ) • ( V X A2{f) ) - Airas A2{f) • ki{r) 



V X y4i(f) ) • ( V X A2{r) ) + 47ras Ai{r) • ^(^ 



(IIL59a) 

E 
111.59b) 



. 



(IIL59c) 

Clearly, these Poisson identities may be obtained also directly from the 
Poisson equations flIII.47aP - flIIL47dp by multiplying through with the appro- 
priate potentials and integrating by parts. But their deduction from the RST 
action principle does better elucidate their meaning for the variational proce- 
dure: obviously, when looking (by trial and error) for those gauge potentials 
which yield stationarity of the action integral VFrst (lII.SGbll . one can restrict 
oneself to those potentials which obey the Poisson identities. It is just with 
reference to this meaning of restrictive conditions that the Poisson identities 
will readily be used in order to set up the RST principle of minimal energy! 

E. Double- Valued Wave Functions 

After the first deficiency (i) of the Ritz-Hartree-Schrodinger approach is 
now eliminated, one can turn to the next critical point, namely the con- 
ventional assumption (ii) that the wave functions must always be unique. 
Indeed we will relax now this presumption and will (as a counterexample) 
admit double- valued wave functions of the type 



(III.60) 
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where {r, {}, 0} are the usual spherical polar coordinates. We will readily see 
that such a more general class of wave functions can generate unconventional 
gauge potentials, namely via the solutions of the Poisson equations fllll.48aP - 
(IIII.48dl) or their non-relativistic approximations resp; and the corresponding 
unusual form of interaction force may then yield energy levels which are closer 
to the conventional Schrodinger predictions than it is possible for the Dirac- 
Fock approach [15]. But clearly, the admitted non- uniqueness (1111. 60|) of the 
wave functions must not imply the non-uniqueness of the physical densities, 
e.g. of four-current fca^ fill. 451) or of energy- momentum density ^^^T^,y flll.50b|) . 
etc. Indeed, it is easy to see that those physical densities are bilinear con- 
structions of and ^ and therefore remain invariant against the change OII1.60p 
of the wave function! 

For a concrete exemplification of those double-valued wave functions, 
one may resort to two basis systems {uJq~^\uJq^} and {cj cuj"'*} of the two- 
dimensional Pauli spinor space which are eigenvectors (with zero eigenvalue) 
of the total angular momentum = Lz + Sz '■ 

JzUJ^^^ = Jz4^ = JA-^^ = 14^ = , (111.61) 

see ref. [B] for the details. Since these basis spinors themselves are already 
double-valued (i.e. ujlj^\r,i!},(f) + 2-11) = —ujl^\r,i!},(j)); etc.), one can decom- 
pose the Pauli spinors ("V±(r) PII-25P with respect to these double- valued 
basis systems as follows (a = 1, 2) 

(")y,+ (f) = (r sm^)-^2 [('^)i?+(r, ^) ■ + (")S+(r, ^) ■ (HI. 62a) 

("V-(^) = -i{rsm^)-^ [('*)^„(r,^9) ■ uj['^^ + (")5_(r,^9) ■ c^i'^j (111.62b) 

(('^V±(r,^9,0 + 27r) = -(")^±(r,^9,0)) , 
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and then both these Pauh spinors and the corresponding Dirac spinors ipai'f^ 
are double- valued in the sense of equation flIIL60|) . provided the wave am- 
plitude ^"'^R±(r,'d) and ^"'^S±{r,'d) are single-valued! Moreover, the latter 
property of uniqueness is transferred also to the Dirac densities ka^{r) = 
{^""^ko (r) , —ka{r)} flll.45p since these appear as bilinear constructions of the 
Pauli spinors ^"'^^±{r), i.e. 

(a) 62 , {a)n2 , (a) 62 , (a)r-2 

7^ f-^ - (III.63a) 
4:7ir sm V 

(III.63b) 

with the azimuthal component '^'^^k^ of the Dirac currents ka{r) being given 
by 

sin^9 f(")^+ ■ ^''^R^ - ('^)^+ ■ ^''^S\ - cos^ (^"^3+ ■ ^"^R^ + ('^)^+ ■ 

(a) L ^ V I V 

27rr sin d 

(III.64) 

Observe here, that through the choice of real-valued wave amplitudes 

, ^"O-t , the radial ('•"^/Cr) and longitudinal {^"''^k^) components of the 
Dirac currents ka{r) do vanish (i.e. '•"^fcr = ^"'''k^ = 0), so that these three- 
currents ka{r) encircle the axis of the spherical polar coordinates {r^'d^cj)). 
Naturally, this symmetry of the three-currents may then be transferred also 
to the vector potentials Aa{r) which thus appear in the following form 

l,(r-) = (")A^(r,^9)-e^, (III.65) 
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from which the magnetic fields Ha{r) can be computed by means of the usual 
curl relation fllll.52bl) . 

Summarizing the properties of the wave functions to be used for the RST 
description of bound states, one first has to mention their double- valuedness 
f lIII.60p and moreover we will assume that the wave amplitudes (")^±,(")^± 
are unique and non-singular (real- valued) functions over space time. But ob- 
serve here that, through this second assumption, the Pauli spinors '^°''^Lp±{f) 
(]III.62aP - (]III.62bp and therefore also the original Dirac spinors V^a(0 (1111.250 
become both singular and double- valued (-w "exotic states"). It should ap- 
pear as a matter of course that such exotic states will imply further un- 
conventional elements of the theory, e.g. the form of the gauge potentials. 
Notice, however, that the observable objects of the theory (i.e. the densi- 
ties of charge, current, energy-momentum etc.) are well-defined and unique 
objects over space-time, the singularities of which (if present at all) do not 
induce any pathological element into the theory. Therefore the wave am- 
phtudes (''^^i, (")5± fllll.62ap - fllll.62bj) as the unique and (mostly) regular 
constituents of the non- unique and singular wave functions ipai'r) flIII.25p 
will appear as the solutions of a well-defined eigenvalue problem. The corre- 
sponding eigenvalue equations are to be deduced from the original eigenvalue 
equations flIII.26aP - flIII.27bp for the double- valued Pauli spinors *^")(^±(r) and 
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appear then in the following form, e.g. for the first particle (a = 1) [7j 

^ + + (^Uo ■ (^)i?. - (^U, fsin^ ■ - cos^ ■ 

or r ov V / 

(III.66a) 

^ - + ^'^^0 ■ ^'^S. + (% fsin^ ■ (1)5+ + cos^ ■ (1)^+) 

or r ov \ / 

h 

(III.66b) 

r or r ov V / 

(III.66c) 

1 a(rW^ + 1^ - (^Uo . - (^U, (sin^ . (i)^_ + cos^ . 
r or r ov V 

n -'"^^ 

(III.66d) 



Mp +Mi 



An analogous set of four eigenvalue equations does apply to the second parti- 
cle (a = 2) which, however, needs not explicitly be reproduced here because 
it can be obtained simply by means of the particle permutation symmetry, 
see ref. [7]. 

F. Unconventional Potentials 

The interesting point with these double-valued wave functions is now 
that they do generate a rather unusual form of the gauge potentials (")y4o(r*) 
and Aa{r) by means of the recipe 0III.48ap - flIII.48dl) . In order to see this 
more clearly, one substitutes the chargg'and current densities ^^^ko{r) 0III.63ap 



and ka{r) (lIII.63bp into those formal solutions of the Poisson equations which 
yields explicitly for the electric potentials in terms of the unique wave am- 
plitudes 

(III.67a) 



4:71 J r'sm-d' \\f—f'\ 



All J r'sin'd' \\f—f'\\ 

(III.67b) 

and similarly for the magnetic potentials '^'^M^(r, -(9). 

In order to estimate qualitatively the new feature of these potentials due 
to the non-sineular wave amplitudes ^''^R±,'-''^S±, it may be sufficient for the 
moment to adopt a simple model ^''^ko^r) for a non-spherically symmetric 
and singular charge distribution being normalized to unity, cf. (III.44I) . i.e. we 
put 

Wfco(r) = ^ = (Anrsin^)-^ . ° gxp -2- (III.68) 

47rr sm v nr^ \ r^J 

with the regular and spherically symmetric charge distribution ^^^ko{r) being 
normalized as follows: 

/ drr (^)fco(r) = - . (III.69) 
Jo ^ 

The corresponding electric potential ^"^^Aq^t) 

.M,fl.a.f-^p^ (111,70) 

J Anr smv \\r — r'\\ 
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will then be found to be also non-splierically symmetric, but it can be 
shown [6] that the binding ability of this potential is supplied mainly by 
its spherically symmetric part (fPly4o(r), say). The latter part may be defined 
by suitable expansion of the denominator ||r — f'\ \ in the integral fllll.TOp . 
see ref. [B] for this method; or otherwise one may substitute the anisotropic 
density ^''^ko{r) (jIII.681) into the RST action principle (lII.56a|) -( ]II.56b|) and 
may then determine the desired isotropic part fPlAo(r) of ^"^^Ao^f) (jIII.701) via 
the solution of the corresponding variational (i.e. Poisson) equation. Resort- 
ing here to the second method it suffices to consider merely the electrostatic 



part Wrsx of the two-particle action integral flll.56b|) with *^^^y4o 
[pUo and = '■^^ko = which yields 



w 



(e) 
RST 



he dr r 



2hc [PUo(r) ■ ^''^koir) 
7r[PUo(r) (^)A;o(r) - 



he 



|VWAo(r)|| 



;ill.71) 



r 

as 



dr 



Thus the electrostatic variational equation (5Wrst = 0) emerges as a spher- 
ically symmetric Poisson equation: 



\ dr"^ r dr J 



o(r) 



vr 



(III.72) 
fUTLGSj) yields 

(III.73) 

This interaction potential, being typical for the exotic quantum states, has 
some peculiar properties: First, it approaches the Coulomb potential (~ — ) 



2 r 

Finally, substituting here the assumed charge density *^'^)A;o('^; 
for the electric potential ^^^AQ{r) 



as A 


'-2^ 




— 1 — exp 




) 


r \ 
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at spatial infinity (r oo) as expected since it is generated by just one 
electric charge unit (i.e. elementary charge). Second, the potential remains 
finite at the origin (r = 0) 

[pUo(O) = ^ (III.74) 



as well as the corresponding electric field strength Ep (IIII.52aP 

c/[pUo(r 



r=0 r? 



(III.75) 



r=o dr 

Moreover, if the length parameter r* tends to zero (r* 0), both the po- 



tential flIII.74p and its field strength flIII.75P approach infinity which says 



that the asymptotic Coulomb form fills then the whole three-space. This is 
clear because, in this limit (r* 0), the charge distribution *-'^-'A;o('") becomes 
pointlike (see fig. 1 of ref. [6j for a sketch of the unconventional potentials). 

However from the physical point of view, the most interesting feature of 
those potentials ^°'^Af){f) due to the exotic quantum states surely refers to 
the fact that they carry a finite energy content, in contrast to the Coulomb 
potential. Indeed, the electrostatic interaction energy E^^^ of both charges is 
given by [5] 

= --^ / d¥||V(Pko(f)||2 , (III.76) 
4vras J 

and when the spherically symmetric approximation [PlAo(r) ( 1III.73I) is sub- 
stituted herein, one finds the following result 



E^^ = -- . (III.77) 



Incidentally, this is just the interaction energy of two point charges sepa- 
rated by the distance which plays the role of a length parameter for our 
model charge distribution ^^^koi^r) flIII.68p . This charge distribution becomes 
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pointlike when the length parameter tends to zero and, clearly, for this 
limit the interaction energy E^^ flIII.77p of both extended charge distribu- 
tions becomes infinite, just as is the case with two point charges of vanishing 
separation (r^, 0). 

Of course, the interaction energy E^^ of the two particles is only a fraction 
of their total energy i?T which must contain also the kinetic form of the 
particle energy. Indeed, this latter form of energy deserves a closer inspection, 
too; and this can be performed most adequately by setting up now the total 
energy functional E^ through adequately exploiting the intrinsic RST logic. 
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IV Energy Functional 



Besides the use of exotic quantum states and their unconventional poten- 
tials, it is necessary to introduce a further new element into the theory in 
order to deal successfully with the energy spectra of the bound systems: This 
refers to the construction of a suitable energy functional i?T, which equips 
the RST field configurations with an energy content being then immedi- 
ately accessible to spectroscopic test. Recall here the fact that, though the 
variational approach flIII.4ip due to the total mass functional Mtc^ flIII.24D 
can be viewed as the relativistic generalization of the Ritz variational prin- 
ciple (lIII.33aP - (]III.33bp . this relativistic approach nevertheless fails to es- 
tablish the gauge field equations and exclusively reproduces the eigenvalue 
equations for the matter fields; see the critical comments (i) and (ii) men- 
tioned above. Therefore it suggests itself to restart from the original notion 
of field energy Ej^ flll.53p - flll.55bl) and to convert this to the wanted energy 
functional E^- 

A. Relativistic Construction 

First, observe that the underlying energy-momentum densities *-°^T^i/ 
and ^'^^T^iy are already specified by equations flII.50ap - OII.52p : and if one sub- 



stitutes therein the stationary form of the matter and gauge fields f llll.laP - 



flIII.2p . one finds the individual energy contributions Ed and Eq flII.55aP 
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(111.55bp appearing in the following form [7]: 

(IV. la) 

= E^' + E^-' . 

(IV. lb) 

Here, the mass equivalents M^'^^c^ and M^i^c^ of the electrostatic interaction 
energy have already been defined previously through equations (lIII.12aP - 
(]III.12bp . Furthermore, the gauge field energy Eq consists exclusively of the 
energy content due to the real gauge field modes A"^, since the complex 
field modes 5^ must be put to zero together with their energy content Eq 
(see the discussion below equation ( III. 261) ). 

Clearly it is very tempting now to consider the total energy functional 
£'t (1II.53|) . with Ey) and Eq being specified by the present equations (llV.laP - 
(jlV.lbp . as the wanted object of our interest. The corresponding variational 
procedure (^-Et = 0) must then be complemented by the former constraints 
of wave function normalization flIII.22p . which had to be applied already in 
connection with the mass functional approach flIII.4ip . By this arrangement, 
one would be led to the following first proposal e!^^ for the desired energy 
functional: 

-^t'' ~ Et + Ad(i) ■ A^D(l) + Ad (2) ■ ^D(2) 

(IV.2) 

= Mtc' + [E^n^ - Mi'^c^ - M^^c^j + E^^^ . 
But observe here that the additional appearance of the electrostatic mass 
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Er 



Et; 



nc 



EAr) ' E^if) + HAr) ' H2(r] 



equivalents and Mfj does spoil the partial success already obtained 

with the mass functional Mtc^ flIIL24p from which the matter eigenvalue 
equations can actually be deduced. The reason is that those mass equiv- 
alents (lIII.12aP - (lIII.13bp do also contain the wave functions ipai'f^, namely 
via the densities ^""^k^^r) flIII.49a|) - (]III.49bp . Therefore we have to eliminate 
again these redundant mass equivalents from our first proposal (lIV.2p . which 
can be achieved by expressing them in terms of the electromagnetic gauge 
field energy E^^^ and E^^^ f HVlbD as follows: 



4^) = M^'^c^ = M^fc^ (lV.3a) 
Eg") = -M/^^^c^ = -M^fc^ , (IV.3b) 



and these relations represent nothing else than the Poisson identities ( 1III.58P 



f lIII.59cp . Thus using this electric coincidence ( llV.SaP in order to eliminate 
the electric mass equivalents from the previous proposal f lTOj) . one 
arrives at the next proposal e'^P: 

eJI) = Mtc' - + . (IV.4) 

This second proposal displays now some pleasant features and therefore must 
be expected to come close to the wanted final result: First, contains the 
matter fields ipai"^ only in form of the mass functional M^c^ and therefore 
the variational equations of the functional must correctly reproduce 
the mass eigenvalue equations (]III.15|) and (1111.17^ ! Second, returning for 



the moment to the matter energy Ej^ (llV.laP and substituting there the 
mass eigenvalues M^c^ (]III.19aP - (lIII.19bp lets the matter energy E-q appear 
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essentially as a sum of single-particle contributions 



2 



Er, = Y,EBia) (IV.5) 

a=l 

with the individual contributions -Ed (a) being given by 

^D(i) = - (M[i]c2 + m/*^)) = Zf,^ ■M^c' + 2Tkin(i) + Mi'^^c^ (IV.6a) 

^D(2) = (M[2]c2 - Mi^^C^) = • M, + 2rkin(2) + M^f^^c^ . (IV.6b) 

This physically plausible result says that the matter energy i?D{a) of either 
particle (a = 1, 2) consists of rest mass energy (first terms) plus kinetic 
energy (second terms) plus magnetic interaction energy (third terms), while 
for these single-particle energies -Ed (a) there appears no electric interaction 
energy. The emergence of the magnetic kind of interaction energy seems to 
be somewhat unreasonable; but this is to be understood as the field theoretic 
counterpart of the minimal substitution {p ^ p ~ ^A) for the conventional 
energy {H = ^) of a point particle moving in a magnetic field H = V x A. 

The third interesting point with that second proposal e!^^^ (1IV.4P is now 
that by use of the explicit form (]III.19al) - (]III.19bp of the mass functionals to- 
gether with the electric and magnetic Poisson identities, cf. ( ]IV.3al) - (11V.3bl) . 
this proposal can be rewritten as the sum of the individual rest mass and ki- 
netic energies plus the gauge field energy of the electric (e) and magnetic (m) 
type: 

E^lll) = ■ Mp + Zf,^ ■ Me C^) + 2 (Tkin(l) + Tkin(2)) + (4'^ - 4"^) • 

(1V.7) 

Here, the validity of the normalization conditions flIII.22p has tacitly 



been assumed and therefore they do not appear explicitly in the present 
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result for e!^^^\ Observe however that the physically reasonable form flIV.7p 
of the energy functional is a consequence of the fact that its preliminary 
form e!^^^ flIV.4p contains the gauge field energies of electric and magnetic 
type with different signs. This important fact is the reason why the double 
counting of the electric term (lIV.3aP in the sum M^c^ (jIII.241) of mass eigen- 
values M[a]C^ (IIII.19ap - (]III.19b|) becomes compensated (see the discussion of 
this effect in ref.s 0ITD]); and then the electric field energy E^^^ appears only 
once in the third proposal E^^^^ flIV.7p . For the magnetic field energy -E^^ 
there occurs an analogous effect since its (negative) double-counting in the 
sum Mtc^ flIV.4p is weakened so that the magnetic field energy E^^ ap- 
pears now in the third proposal (]IV.7p with the opposite sign relative to 
its electric counterpart -E^'*! This circumstance however does not influence 
the lowest-order approximation of the atomic energy levels because these are 
dominated by the electric interactions. Nevertheless for the higher-order ap- 
proximations, the negative sign of the magnetic term will leave its imprint 
upon the predictions and therefore must receive confirmation or rejection by 
the observational data (see below). 

The final step for the construction of the wanted energy functional must 
now be based upon the somewhat amazing circumstance that the third 
proposal (1IV.7P would numerically produce the same energy upon an 
exact solution of the RST eigenvalue problem as does the original func- 
tional Et flll.53p - flll.55bp . too. The reason is that the transcription of Et 
to i?,^"'* relies exclusively upon the use of the Poisson identities which, how- 
ever, do automatically hold for any exact solution of the RST eigenvalue 
problem! Nevertheless, this third form i?-^"^ f]IV.7p cannot be used for the 
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deduction of the mass eigenvalue and Poisson equations as the corresponding 
variational equations since this functional fllV.Tp contains no coupling at all 
between the matter fields ipa{r) and the gauge fields ^"'^Ao{r) , Aa{r) . Indeed, 
the coupling of matter and gauge fields has been eliminated on the way from 
the orig inal Et ( 111531) to the present E^!"^ fUVTl) via the Poisson identities, 
albeit under simultaneous preservation of the numerical value of the energy 
functional. 

Therefore it finally becomes necessary to restore that lost coupling of 
matter and gauge fields, again under preservation of the numerical value 
of the energy functional. Naturally one expects that such a restoration of 
the desired coupling must be performed with regard again of the Poisson 
identities which thus have to take over the role of constraints for the varia- 
tional procedure (see also ref. [6]). In this sense, one resorts to the method 
of Lagrangean multipliers with respect to both the wave function normal- 
izations (1III.22P and the Poisson identities (1III.58p - (lIII.59cP : and thus one 
complements the third proposal E^^^^ (\lV.7\i to the final result E^ in the 
following way: 

E^ = + f: (Ad(„) ■ N^ia, + Acg ■ iVog + Ao[:)^ ■ No^^^) . (IV.8) 

a=l 

Here it is now a standard exercise to convince oneself of the fact that the 
variational equations due to this functional E^ actually are just the mass 
eigenvalue equations flIII.151) and flIII.171) together with the electric and mag- 
netic Poisson equations flIII.47a|) - flIII.47dl) . provided the Lagrangean matter 
multipliers Ad (a) are given in terms of the mass eigenvalues as shown 
by equations (IIII.23aP - flIII.23bp and furthermore the gauge field multipli- 
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ers Ag(q™'' are specified as follows (a = 1, 2) 

AoS? - ^ . (.V.9b) 

Thus collecting all the partial results, the ultimate form of the wanted energy 
functional is the following: 

Et = Zfl) ■MpC^ + ■ Me + 2 (Tkin(l) + Tkin(2)) + (4'^ " 4"^) 

^ a=l 

(IV. 10) 

The practical usefulness of this ultimate energy functional E^ ( IIV.IOP 
refers mainly to those situations where the RST eigenvalue problem cannot 
be solved exactly so that one is forced to look for approximate solutions 
(which will be mostly the case). But fortunately, a convenient approxima- 
tion method is now at hand in form of the energy functional -Et, so that one 
can test certain trial functions for the Dirac spinors ipai'f^ and for the gauge 
potentials ^"'^Aq^t) , Aa{f) . These trial functions will depend upon some pa- 
rameters {bk) so that, after substitution of the trial functions into the energy 
functional E't (1IV.10|) . the latter becomes an ordinary function of the ansatz 



parameters bk'- -Et = ET^{bk)- Finally, looking for the minimally possible 
value of that function E-^{bk) yields a more or less good approximation for 
the wanted energy eigenvalue of the RST eigenvalue problem. For an exam- 
ple of this type see ref. [6]. However a further improvement of this general 
approximation procedure may be achieved by not trying some independent 
functions for the gauge potentials '^"^y4o(r) and Aa{r) but by trying merely 
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for the (normalized) wave functions V'a(r) then calculating (exactly) the 
associated gauge potentials from their Poisson equations flIII.47aP - flIII.47dp . 



preferably in form of the solutions fllll.48al) - fllll.48dl) . Clearly, through such 
a procedure the Poisson identities flIII.58p - fllll.59cp will be satisfied exactly, 
though the associated solution of the RST eigenvalue problem is an approx- 
imation. But the advantage is here that all the constraints (second line on 
the right of equation fjlV.lOp ) can be omitted and one can concentrate upon 
the physical terms (first line) which effectively is i?,^"^ f llV. 7p . Thus it will 
be sufficient to look for the minimum of the corresponding function E!^^^\bk)- 
Subsequently we will exemplify this procedure by means of the positronium 
groundstate. 

As a preparation of this groundstate treatment, it is very instructive 
and convenient to specify the functional E'-^"'* in terms of the wave ampli- 
tudes ^"^Ri and ^"^S^ fllll.62ap - fllll.62bp . First, the mass renormalization 



factors Z'^^s^ fllll.lip are found to be of the following form 

Zf,) = i y d'f(^^-^Rl + ^-^Sl - ^''^Rl - ('^^S^) , (IV.ll) 

where the unique wave amplitudes ^°'^R±, ^"'^S± are assumed (for the sake of 
simplicity) to depend only upon the radial (r) and longitudinal {i)) vari- 
ables: (")^±(r,^9),(")^±(r,'^); and the remaining two-dimensional volume el- 
ement (i^r is then given in terms of these variables as 

(fr = rdrd'd. (IV. 12) 

Next, the kinetic energies Tkin(a) f HIl.20aP - fllll.20bp of both particles 
(a = 1, 2) are found to split up into the radial (T^) and longitudinal (T^) 
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part, i.e. 

^kin(a) = Tr(a) + T^(a) (IV. 13) 

with the radial part being given by 



Trio) -{-I) -ljdr\ R.-^^ ^ 



dr r dr 



and analogously the longitudinal part by 



d'd 



(IV. 14) 



(IV.15) 



Furthermore, the electric and magnetic field energies e'I^ and E*^^ firvlb|) 



read in terms of the gauge potentials ("Mo(^) and '^"^^(r) ( 1111.651) 



^(e) he f,,^ . J d^'Kjr) g(%o(f) , 1 g(%o(f) 9(%(r) 
Et, = / a r r sm w ■ h 



R 2a„ / \ dr dr d'd d'd 



^ 2as / '^"^ Kdr"-' dr 



+ 

sm 



(IV. 16a) 
(IV. 16b) 



Finally, both kinds of constraints, i.e. the normalization conditions (111.440 
and the Poisson identities (IIII. 58^ - 0111. 59cl) . must also be rewritten in terms 



of the wave amplitudes; but it is not necessary to reproduce this here because 
for the subsequent treatment of the positronium groundstate we will use trial 
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functions satisfying a priori all those constraints and therefore we can rely 
directly upon the truncated functional i^r^""* fllV.Tp without loss of accuracy. 



Now in order to support the confidence in the established functional 
-Et fllV.lOp . one can look for both the mass eigenvalue and Poisson equations 
in terms of the wave amplitudes ^"'^R±, by carrying out the variational 
procedure {6Et = 0) just with respect to these wave amplitudes and gauge 
potentials. Clearly, one will then actually recover the former mass eigenvalue 
equations (lIII.66aP - (]III.66d|) together with the Poisson equations (lIII.47aP - 



f lIII.47dD . For their magnetic part flIII.47cD - flIII.47dD one may resort for the 



moment to the special case of circular flow around the z-axis, cf. flIII.65P : 



and in this special case the magnetic Poisson equations for the azimuthal 
component ^""^A^ read 

sm v 

A(^M^ - = Ana^ , (IV.17b) 

sm v 



with the circular current components '■"-'fc^ being given by equation flIII.64p . 

B. Non-Relativistic Approximation 

For a first practical test of the present construction of an RST energy 
functional fllV.lOp . it may be sufficient to restrict oneself to the non- 
relativistic approximation. Clearly, if such an approximation would fail to 
meet with the well-known results of ordinary non-relativistic quantum me- 
chanics, one would not try to further elaborate the corresponding relativis- 
tic situation. Fortunately, the subsequent demonstration by means of the 
positronium groundstate points just into the other direction: The conven- 
tional groundstate energy can be exactly reproduced by an appropriate trial 
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function for the non-relativistic limit of the RST functional (IIV-IOP - 
or E^"^ fHv:7|) . resp. In order to find the desired non-relativistic limit of the 
functional e!^^^\ it is merely necessary to look for the non-relativistic forms 
of its constituents, i.e. rest mass and kinetic energy and the field energy of 
the electric (E^^) and magnetic kind (E^^). 

Naturally, the non-relativistic situation becomes even further simplified if 
one restricts oneself to the spherically symmetric approximation by neglecting 
the magnetic interactions. As a matter of course, the electric fields Ea{r) can 
easily be visualized to be spherically symmetric ( hedgehog configuration) 
in contrast to the magnetic fields which mostly obey a dipole (or higher) 
symmetry. Therefore it is favorable to start with the spherically symmetric 
configurations of the purely electric type. 

Turning here first to the mass eigenvalue equations (]III.66al) - (]III.66dp . 
one usually assumes that the "negative" Pauli components ^'^-'i?. , are 
much smaller than their "positive" counterparts and ^'^o^, so that 

the non-relativistic form of the mass eigenvalue equations is obtained by 
simply eliminating those negative components ^^^R-,^'^^S- [7\. The residual 
eigenvalue equations for the positive components of the first particle (a = 1) 
do appear then in the following form: 



2M„ 




dr 



1 9^(1)/?, 




+ 



1 d-^ws. 



+ nc(%(f) ■ (^)A 



^1 



8(1) 



(IV. 18a) 



+ /ic(%(r-) ■ (1)5+ = Es(i) ■ (')^- 



(IV. 18b) 
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Here, the magnetic interactions are neglected ('•"M^ — > 0) together with the 
relativistic effects because both phenomena are mostly of the same (small) or- 
der of magnitude. Moreover, the (conventional) non-relativistic Schrodinger 
eigenvalues -Es(a) are defined as in equations flIII.30aP - flIII.30bl) . The case 



of the second particle (a = 2) is not written down because it can easily 
be supplied by means of the particle permutation symmetry (1 ^ 2), see 
ref. [7]. However the important points with the non-relativistic eigenvalue 
equations flIV.18aP - fjIV.18bp refer now to the facts that (i) they are not of the 
usual Schrodinger form flIII.32aP - flIII.32b|) and (ii) the spin-up (~ R+) and 



spin-down (~ 5*+) configurations are decoupled. The latter circumstance ad- 
mits us to conceive either of the two single-particle spins to point definitely 
into the positive or negative z-direction and their combination to the para- 
and ortho-states of the two-particle system will then intuitively be evident. 

Naturally, one expects that these non-relativistic eigenvalue equations, 
such as f lIV.18ap - fjIV.18bp . should emerge as the variational equations due 
to the non-relativistic approximation {e!^\ say) of the original RST energy 
functional f llV.lOp . Indeed, one is easily convinced that this supposition 
is true; namely the elimination of the negative Pauli components ^"•'i?^, ^"o. 
from the relativistic kinetic energies Tkm{a) flIV.13p - flIV.15p yields [6J 

(Zfi ) - l)Mp + 2Tkin(i) ^kin(i) + ^w(i) (IV. 19a) 

{Zf^^ - l)Me + 2Tkin(2) ^kin(2) + ^W(2) (IV.19b) 
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with the non-relativistic kinetic energies i?kin(a) being given by 



kin(l) 




kin(2) 



and the "winding energies" -Ew(a) being given by 



(IV.20a) 
(IV.20b) 
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W(l) 
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4M. 
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d^'^R^ 9(1)^4- 9(1)5-. d^'^R, 



dr di!) 



9r 



dr 



dr 



(IV.21a) 
(IV.21b) 



Observe here that, for the kinetic energies i?km{a)) we made use of the non- 
relativistic decouphng of the spin-up and spin-down components and thus 
adopted the first spin (llV.lSap pointing in the positive z-direction {-^ ^^^R+) 
and the second spin fllV.lSbl) in the negative z-direction (-w (2)5+). Clearly, 
the other combinations of the spin directions {(i)i?+, (2)/?+}, {(i)5+, (2)-R+}, 
{(1)5+, (2)5+} are equally well possible, see below for the para- and ortho- 
configurations. Fortunately, through this choice of definite spin directions for 
any particle, the winding energies (lIV.21aP - (]IV.21bp become zero so that one 
can restrict oneself to the kinetic energies (]IV.20a|) - (]IV.20b|) alone. Of course 
the rest mass energies need not be taken into account for a non-relativistic 
treatment and therefore have been omitted, cf flIV.19al) - flIV.19b|) . 

Next, the non-relativistic form of the electric and magnetic gauge field 
energies and e'^^^ remains the same as in the relativistic case, i.e. in 
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terms of the static gauge potentials ('^Uo(r^ and Aa{r) (cf. (lIII.52aP - (lIII.52bp : 



4^) = j d'r (V(%(f)) ■ (V(%(f)) (IV.22a) 

= j d'r (V X Ai) • (V X A2) . (IV.22b) 

But for the exphcit calculation of the non-relativistic potentials ^"•^AQ{r) 
and ^"^0(7^) from the Poisson equations flIII.47aP - flIII.47bD and fllll.47cp - 



flIII.47dp one will use the corresponding non-relativistic approximations for 
the charge and current densities, cf. f lIII.63aP and f ]III.64p 
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(lV.23a) 

47rr sm v Aur sm w 
l-nr sm w 27rr 

Observe here again that, in the contrast to the charge densities *^'^-*fco (r), the 
current densities '^'^\^{r) are built up by both the positive {^"'^R+, ^""^S^) and 
negative (^"■'-R_, ^""^S^) wave amplitudes while, properly speaking, the nega- 
tive amplitudes ^""^R^, ^""^S^ should be neglected against their positive coun- 
terparts ^^^R+, ^^^S+ for the non-relativistic limit. This demonstrates that it 
appears somewhat inconsequent to retain the magnetic (i.e. spin-spin) in- 
teractions for the non-relativistic approximation because their order of mag- 
nitude may be the same as that of the other dominant relativistic effects. 
Nevertheless we will not drop the magnetic effects for our non-relativistic 
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approach because one can still deal with the effect of ortho-para splitting of 
the energy levels from a more qualitative viewpoint. 

In this sense one has to renounce on the inclusion of the magnetic ef- 
fects for the purpose of deducing the non-relativistic eigenvalue equations 
f lIV.18ap - flIV.18b|) from the desired non-relativistic version E^?^ of the orig- 
inal functional Et ( 1IV.10|) . Consequently one drops also the magnetic con- 



straints Ng^/^\ ( 1III.59bl) - (]III.59cl) from the latter functional and retains only 



the electric constraints ^G(a) (lIII-58p - (lIII.59ap which do then appear in the 
following form: 



(e) 
(1) 



AT 



^G(2) 



+ a. 



(%o(r) ■ (2)^^ 



r sint? 



(IV.24a) 



(IV.24b) 

where the non-relativistic approximations (lIV.23aP - (]IV.23bp of the charge 
densities ^""^ko (r) have already been respected. Clearly, the latter approx- 
imations for the charge densities must also be used for the constraints of 
wave function normalization flIII.22p which then appear in their following 



non-relativistic forms Nj) 



(0). 
(a)- 



AT, 



D(l) 



(2) 



. 



(IV.25a) 
(IV.25b) 



Finally, collecting all the non-relativistic approximations and applying 
again the method of Lagrangean multipliers lets appear the wanted non- 
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relativistic approximation e!^^ of the original functional i?T fllV.lOj) in the 
following form: 



E. 



(0) 
T 



-E'kin(l) + E] 



kin (2) 



+ E 



(e) 



2 2 (IV.26) 

a=l a=l 

Here it is again a nice consistency check to convince oneself of the fact that the 
usual variational procedure {6e!^^ = 0) actually does reproduce the claimed 
non-relativistic forms fllV.18ap - (IIV.18bp and (llll.47ap -( iril.47bp of the mass 
eigenvalue and Poisson equations. The non-relativistic multipliers As (a) turn 
out as the conventional Schrodinger energie eigenvalues 



A 



S(a) 



-E> 



S(a) ) 



(IV.27) 



which compares to the analogous result (1III.40P of the Ritz-Hardy-Schrodinger 
approach. Clearly according to our present choice of the negative z-direction 
for the second particle spin, the second eigenvalue equation (jIV.lSbp for the 
first particle must be replaced for the present situation by 



2Me 



1 d 



+ 



1 92(2)^4 



he (^Un ■ (2)^4 



"-^8(2) ■ ^^^'S'^ 



(IV.28) 

for the second particle. Furthermore, the non-relativistic approximations of 
the electric Poisson equations do now appear as the following variational 
equations: 



a. 



r sin'i? 
r siwd 



(IV.29a) 
(IV.29b) 
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The important point here is that, despite the many similarities between 
the non-relativistic hmit of RST and the conventional Ritz-Hartree-Schrodinger 
approach, there are also characteristic differences of both approaches which 
are in favour of RST. This will readily be demonstrated by considering a 
numerical example. The main difference refers to the gauge potentials, e.g. 
those of the electric type (lIII.48al) - (]III.48bp . which by means of the non- 
relativistic approximations pV.23aP - (lIV.23bp for the charge densities ^""^ko (r) 



appear as 

.^Mrl^"ff^-^p^ (.V.30a) 

Air J r sm v' \\r — r'\\ 

^-^/^.^f<^. (.V.30b) 

Arc J r smij' \\r — r'\\ 

Obviously, these gauge potentials due to the exotic states must be more 

singular as the Hartree potentials flIII.50ap - flIII.50bp which are due to the 



non-singular Hartree wave functions fai'f^) f1III.3ip . see the example flIII.73p . 

C. Magnetic Interactions 

In the contrast to the electric fields Ea{r) flIII.52ap . the magnetic fields 
Ha{r) flIII.52bp cannot obey the SO (3) symmetry because they have dipole 



character rather than monopole character like their electric counterparts. 
Therefore it will become necessary to apply more complicated approximation 
techniques; but fortunately it is not necessary to explicitly solve the magnetic 
Poisson equations flIII.47cP - fjIII.47dp for the three- vector potentials Aa{r), 



e.g. in form of the special solutions flIII.48cP - flIII.48dp . Rather it is sufficient 
to determine the magnetic fields Ha{r) {a = 1,2) directly from the Abelian 
Maxwell equations 

Vx#, = 47raJ„, (IV.31) 
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which is the three- vector form of the relativistic versions (III.27ap - flII.27bp . 
Nevertheless, one has to insist on the existence of the corresponding vector 
potentials Aa{f) fllll.52bl) . namely in order that the magnetic Poisson identi- 
ties fllll.59bp - pil.59cp can be satisfied and thus the corresponding magnetic 
constraints in the energy functional E-r fllV.lOp can be dropped. Indeed, in 
the latter case one can restrict oneself to the physical terms of the energy i^x 
(i.e. the first line on the right-hand side of fllV.lOP ). where the magnetic inter- 
action energy is then simply given by E^^^ fllV.lbp in terms of the magnetic 
fields Ha themselves. But clearly if there is no difficulty with the determi- 
nation of the vector potentials Aa{f) fllll.52bp directly from their Poisson 
equations, one may calculate the magnetic interaction energy ij") also m 
terms of these vector potentials Aa{r) as shown by equation fllV.22bp . 

Following here the first path (i.e. determination of the magnetic fields Ha{r) 
directly from the Maxwell equations f]lV.3ip with omission of the vector po- 
tentials Aa{r)), one additionally has to impose the conditions of vanishing 
sources 

V • Ha{r) = , (IV.32) 

in order to ensure the existence of the vector potentials Aa{r). Since the 
Maxwell currents ja{^) are connected to the Dirac currents ka{r) by equa- 
tions fill. 43aP - fill. 43bp . the Abelian Maxwell equations fllV.3ip read in the 
component form of the spherical polar coordinates {a = 1,2) 



dr ' d-d 



47ra,(")A;^, (1V.33) 



where the azimuthal components ^"'^k^ of the Dirac currents ka{j) are speci- 
fied by equations f llV.23cP -f llV.23dp . However for the present non- relativistic 
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limit, the "negative" Pauli wave amplitudes ^°'^R^,^"'^S^ must be approxi- 
mately traced back to their "positive" counterparts '•"^-R+, ^""^8+ through [HI [7] 




2Mp c (9r 

^ 2MpCy dr r d'd j 2M^c dr ^ ■ > 

and analogously for the second particle (1 — ^ 2; Mp — > Me ). This approxima- 
tive procedure recasts the non-relativistic current components '•"-'A;^ pV.23cP - 
( 1IV.23dl) to the following form: 

Wk. ~ (IV.35a) 

* 2Mp 27rr ^ ^ 

^ 2Me 27rr ^ ^ 

Evidently, the currents ka{r) become singular at the origin (r — > 0) unless 
the wave amplitudes *^^'*-R+, ^^^5+ or their derivatives do vanish for r — > 0. 

Of course, these singular currents will imply a corresponding singular 
behavior of the magnetic components ^"''^H^ and ^""^H^ as solutions of the 
Maxwell equations fllV.SSp . In order to elaborate this singular behavior in 
some more detail, one splits off the short-range magnetic field H'^{f) from its 
far-range counterpart which can always be taken as a gradient field, i.e. we 
put 

H,{r) = H'^ir) + Vr^(„) . (1V.36) 

The inclusion of such a gradient field is necessary in order to satisfy the 
divergence relation (1IV.32P which yields for the magnetic potentials rj(^a) 

Ar^(,) = -(v-H'\ . (IV.37) 
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Thus the magnetic Maxwell equations (lIV.3ip do fix only the curl of the short- 
range fields H'^if) and leave the determination of the magnetic potentials ri{a) 
to the Poisson equation flIV.371) . 

Observe also that the superposition of a magnetic potential to the 
original magnetic field H'^ij-) flIV.36p . so that the resulting field Haif) gets 
vanishing divergence (lIV.32p . does not only ensure the existence of a vector 
potential Aa{r) for Ha{= V x Aa) but additionally implies a further effect 
which refers to the magnetic field energy 

^Jr^ fUVTb]) . Indeed the latter 
object reads by use of the superposition flIV.36|) 

= j d^r {h[ 'H'^ + H[' V77(2) + ^2 • V77(i) + V77(i) ' ^V(2)) ■ 

(IV.38) 

Here it is easy to see that for given fields H'^{r) this functional of the magnetic 
potentials rj(^a) is stationary just upon the solutions of the Poisson equations 
()IV.37p ! Thus the introduction of the magnetic potentials //(a) does not only 
guarantee the existence of the vector potentials Aaif) but it additionally lets 
the magnetic energy functional appear stationary (independent of the 
stationarity of the total functional -Et)- 

The ansatz for the components of the short-range fields H'^ is now (a = 1, 2): 

(")ff;(r, ^) = (")if;(r, ^) = (IV.39a) 
(»)i/;(r,^9) = /i„(r) , (IV.39b) 

i.e. the integral lines of the short-range fields H'^ = ha{r)e^ are circles 
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r = const, in the two-planes = const, with the center located at the ori- 
gin r = 0. Clearly, such a field must necessarily be singular along the z-axis 
but this singularity does not contribute to the magnetic field energy 

The magnetic ansatz functions ha{r) are linked to the azimuthal compo- 
nents (")fc<^(r) fllV.35ap - flIV.35bD via the Maxwell equations fllV.331) . which 
do apply also to the short-range components ^""^H'^, ^'^^H!g, yielding 

i A (rhi) = Ana, ^^^ks (IV.40a) 
r dr 

-:^(r/i2) = -47ra,(2)^ (IV.40b) 
r ar 

But since the current components ^'^^k^ are just of that specific form flIV.35aP - 
flIV.35bp . the solut ions of the present differential equations f lIV.40aP - flIV.40bp 
are easily found as 



2 



hi{r) = ^ (IV.41a) 

2Mp c r 



ash 



i2 



i.e. the short-range fields H'^{r) = {^"^ H'^,^"^ H'^} fllV.39al) - flIV.39bD can be 
directly traced back to the non-relativistic wave amplitudes ^^^R+, ^^^S+. 

Unfortunately, the determination of the long-range magnetic potentials 77(a) 
from their Poisson equations flIV.37p is technically somewhat more compli- 



cated. It is true, the source of the short-range fields H'^{r) is relatively simple 

— * — * P0I" 7/ 

W-K{r) = K{r) . (IV.42) 

r 

Thus the corresponding standard solutions for the magnetic potentials //(a) 
are given by 

d?r' cot'd' ■ hair') 



Via) ^ 



, 1 f d-'r' cot'd' -hair') , 

47r y r I |r — r 'I I 
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I ryl'2i 2 V I ryf'2i 



(IV.44) 



where the radial functions hair) are specified by equations (lIV.41aP - (lIV.41bp . 
However it seems here that the integral cannot be calculated in terms of 
analytic functions, not even for the simple exponential trial form fllll.GSP 
for the wave amplitudes '•^■'5'^. Therefore one will be forced to apply 

more or less effective approximation methods. For the present purpose one 
expands the denominator in the integral (1IV.43P as follows 

1 _ 1 

Here the first (i.e. monopole) term does not contribute to the magnetic po- 
tential rji^a) ( 1IV.43I) so that we may be satisfied in the lowest order with the 
dipole approximation (second term). Thus the magnetic potential rjf^a) be- 
comes then in this dipole approximation 

^ (^)r/,(r, d) = -r cos ^ / dr' r'^ " 3 • (IV.45) 

The dipole character of this result becomes evident from its asymptotic 
behavior (r 00): 

(^)r7a(r 00, t9) = - ■ ^ r dr' r'^ha{r') . (IV.46) 
4 Jo 

If the preceding results ( 1IV.41al) -( lIV.41b1) are used here, with observation of 

the non-relativistic normalization conditions (1III.69I) reading explicitly, e.g., 
for the second particle 

2 

dr r ^^^Slir) = - , (IV.47) 

71 

then the second magnetic potential flIV.43P appears in the asymptotic re- 
gion (r 00) as 

'°W-oo,.) = ±.f (.V.48) 
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Apart from the dimensional factor (^), which is due to our use of ge- 
ometric units for the potentials and field strengths (see equation flll.ll|) ). 
the present result flIV.481) for the asymptotic magnetic potential is the usual 
one for a magnetic dipole which however carries only half of a Bohr mag- 
neton /iB(=i= 2Wc^- Clearly, this is a further unconventional feature of the 
exotic wave functions V'a(^); namely besides their singular character for r — 
( ]III.62al) - (]III.62bp . their doubled- valuedness (lIII.GOp . and their integer spin 
eigenvalue (lIII.6ip . 

It must be stressed, however, that these exotic states do not induce any 
pathological feature into the theory, neither with respect to the electric field 
nor for its magnetic counterpart. As a brief demonstration one may inspect 
the magnetic fields in the vicinity of the origin (r — » 0). First, the scalar 
magnetic potential 7^(2) is rewritten as 

(^)r^2(r,^) = |^-^^?2(r), (IV.49) 

with the dipole screening factor g2{r) being given by 

g2[r) = vr -r / ar - 3 — > 1 , (IV.50) 



as h Jq ^/W^ 



12- 



so that the asymptotic dipole behavior (]IV.48P is immediately manifest. (For 
the first particle, a = 1, the same arguments do hold with merely the electron 
mass Me being replaced by the mass Mp of the positive particle). But the 
crucial point with the scalar magnetic potentials ^^^r]2{r-,'&) flIV.45P is now 
that they do not diverge at the origin. One is easily convinced of this assertion 
by tentatively substituting for the magnetic ansatz function h2(r^ its non- 
relativistic form (lIV.41b|l with the wave amplitude ^'^'Sj^ being deduced from 
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the charge density '^*^A;o(^) (1III.68P as 



vrr? ' \ r* 



(2)^+(r) = J^-exp(--) . (IV.51) 



This then yields for the dipole screening factor (72 (r) fllV.SOP 



/ 2 \ 2 roo r exp I - 

,.3.y^ "r'^^, (IV.52) 

which by substitution of the integration variable r' 

p^- (IV.53) 
r 

adopts the following form 

/2\2 /-oo p.exp(-|^-p) 

However in this form, it is easy to see that in the vicinity of the origin (r 0) 
the screening factor (72 (r) looks as follows 

^72(r)^ f-V-r^ (IV.55) 



J*. 

and thus yields a finite value of the magnetic dipole potential ^^^r]2{r,'d) 
pV.49p around the origin. This compares to the analogous behavior of the 
electric potential ^^^Ao{r), see the discussion below equation flIII.73p . 

But once it is guaranteed that the vector potentials Aa{r) do really ex- 
ist, one can use this fact in order to recast the magnetic interaction en- 
ergy f llV.lbp in a new form which exclusively is based upon the Dirac 



currents ka{r): 
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Indeed in order to arrive at this result, one merely has to substitute the 
vector potentials Aa{r) in the magnetic mass equivalents Mj™^c^ or Mj™''c^ 
of the Poisson identities flIV.3bl) by the formal solution flIII.48cl) - flIII.48dp . 
Furthermore, both Dirac currents ka{r) are of the azimuthal form flIII.63bp 
with ^""^k^ being specified by equations flIV.23cp - flIV.23dl) . which in their non- 
relativistic form appear as shown by equations f]IV.35a|) - f]IV.35b|) . Thus the 
magnetic energy does finally emerge in the following form for identical 
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(IV.57) 

Recalling here the fact that for the groundstate both particles must be in the 
same quantum state (apart from the spin direction), one puts 



and thus the magnetic energy flIV.57p becomes 
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(IV.59) 



Here it is reasonable to assume that the optimal value of the variational 
parameter will be found of the order of magnitude of the Bohr radius ob ; 
and this implies that the magnetic interaction energy E^^ is smaller than its 
electric counterpart E^^ flIII.77p by the factor cts ^ ^ 10"'^. This is also the 
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order of magnitude of the other relativistic effects; and therefore the ground- 
state energy difference of ortho- and para-positronium cannot be expected 
to be properly predicted by the present purely magnetic result flIV.59|) (see 
below) . 
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V Positronium Groundstate 



For the situation where both particles masses are identical (Me = Mp = 
M), it is reasonable to assume that both the first (positively charged) par- 
ticle (a = 1) and the second (negatively charged) particle (a = 2) do al- 
ways occupy physically equivalent states. According to this assumption, 
the positronium energy spectrum is expected to be essentially a one-particle 
spectrum which is in perfect agreement with the observational data [S]. The 
conventional classification of the positronium energy levels relies on the com- 
position law for angular momenta so that the groundstate appears as the 
doublet I^Sq and l^^i corresponding to whether the total spin S is zero 
{S = Si — S2 = 0; para-positronium) or is unity {S = Si + S2 = 1; ortho- 
positronium) , see e.g. ref. [16]. However in RST as a fluid-dynamical theory, 
it is more adequate to base the classification upon the relative orientation of 
both magnetic fields Ha{r) rather than upon the angular momentum compo- 
sition law which is adequate for the conventional tensor product of Hilbert 
spaces but not for the present Whitney sum of single-particle bundles. 

But in any case, the inclusion of the magnetic (i.e. spin-spin) interac- 
tions is an additional complication; and it is therefore convenient to first 
simplify the problem by neglecting the magnetic interactions completely and 
considering the residual problem alone (-^ "electrostatic approximation"). 
In the non-relativistic conventional theory, this truncated problem is then 
described by the two-particle Hamiltonian Ho, (]III.34|) and can be solved ex- 
actly by introducing the relative and center-of-mass coordinates [17]. The 
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corresponding conventional groundstate energy -Eolcon is then easily found as 

^olcon = -7— = -l^s ^-Mc^^ -6, 80 [eV] (V.l) 
4 Ob 4 

where ob (= h'^/Me'^) is the Bohr radius and as (= /he) is the fine structure 
constant. 

Indeed, this result (IV. ip is nothing else than the conventional hydrogen 



groundstate energy due to a fixed nucleus, with merely the electron mass M 
being replaced by the reduced mass M/2 due to the comoving positron. Natu- 
rally, one will demand from any new theory of quantum matter that it should 
reproduce this standard result (IV. ip in its lowest order of approximation; and 
afterwards one may proceed to compare the higher-order predictions of the 
various theoretical approaches. Therefore we will now first clarify the way 
in which the standard result fIV.ip emerges in RST, and afterwards one can 



turn to the magnetic effects as small corrections of the electrostatic results. 
Amazingly enough, we will recover just the standard result (]V.1|) as an ap- 
proximate (i.e. variational) solution within the RST framework, namely by 
resorting to the RST principle of minimal energy 

A. Electrostatic Approximation 

Reasonably, the electric properties of both particles may be adopted to 
be approximately independent of the different magnetic arrangements. Thus 
the Dirac densities '•"•'/co (r) can be assumed to be the same for both particles: 

(i)A;o(f) = ^%o{r) = ^% (f) . (V.2) 

Furthermore, since these charge densities generate the electric potentials 
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^"^y4o(r) according to the Poisson equations (lIII.47all - (]III.47b[) . both poten- 
tials can differ at most in sign, i.e. 

(%(r) ^ -(%(rl = ^^%{r) , (V.3) 

with the common potential ^^'^Aq^t) obeying the Poisson equation 

A(PUo(f) = -Ana, (r) . (V.4) 

If we resort here to the non-relativistic approximations ( lIV.23al) -( lIV.23b|) and 
tentatively put for the non-relativistic wave amplitudes, cf. f llll.GSp . 



(2)^+ = (1)^+ ^ R{r) = ^(%o (r) = y ^ ■ exp (^-^^ (V.5) 

then we just recover the former model potential [Ply4o(r) ( 1111.731) as the com- 
mon potential (IV. 30 for both particles. And correspondingly, their electro- 
static interaction energy E^^^ (llV.22aD is then just given by equation (1111. 77^ . 



Here it is important to remark that, due to the non-relativistic trial 
function R(r) (IV. Sp . the potential tPlylo(r) is the exact solution of the non- 
relativistic Poisson equation (1111.721) and therefore the Poisson constraints 
( llll.58P - (llll.59al) are exactly satisfied in their non-relativistic form (llV.24aP - 



( llV.24bp . Furthermore, it is easy to see that the non-relativistic trial func- 
tion R{r) (IV.SP actually obeys the normalization conditions (lIV.25aP - (llV.25b[) . 
Thus both constraints for the (second line of the) non-relativistic energy func- 
tional e!^^ (11V.26P are automatically satisfied by our spherically symmetric 
trial wave amphtude R{r) (IV.SP : and therefore one is concerned solely with 
the physical contributions (first line) to the energy functional. 

However, since the field energy E^^ (as the interaction energy of both par- 
ticles) is already specified by equation (1111. 77p . one is left with the determi- 
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nation of the non-relativistic kinetic energies i?kin(a) (HV-20aP - (lIV.20b|) . Ob 



serving here the spherical symmetry of our trial function fIV.SP together with 



the fact that both kinetic energies must be identical (i.e. -Ekm(i) = E\^m{2)), 
one arrives at the total kinetic energy Ey^i^ as 

h"^ C f dR(r) \ ^ h"^ 

£u,„ = + £k,„,, = ^ J d'r (^ j = ^ . (V.6) 

Consequently, the value of the non-relativistic functional eI^'^ upon our spher- 
ically symmetric trial function R{r) (IV. 51) becomes the following ordinary 
function E^\r^) of the ansatz parameter r^,: 



According to the established principle of minimal energy, the positronium 
groundstate energy Eq in the spherically symmetric approximation is given 
by the minimal value of this function E!^\r^) (IV. 7p . i.e. 



2 



(0) 



Eq = E^ 



e2 



and this minimum occurs for the value r^, =^ 2aB of the ansatz parameter r^,. 
Thus the non-relativistic approximation of the RST principle of minimal en- 
ergy yields just the conventional Schrodinger value (IV.ip for the positronium 
groundstate! This, however, is surely an amazing result in a two-fold way: 

(i) Despite the very different mathematical structure of both approaches 
(Whitney sum vs. tensor product) the present RST prediction (IV. Sp 
coincides exactly with the conventional Schrodinger prediction ( IV. ip . 



(ii) However, in contrast to the conventional prediction, which is adopted 
to be exact within the standard framework of quantum mechanics, the 



corresponding RST prediction (IV. Sp is based upon the choice of an 
appropriate trial function, cf. flV.51) . and therefore is an approximate 
result within the RST framework. Thus the interesting question arises 
how close the exact (but non-relativistic) RST prediction would come 
to the exact conventional prediction fIV.ip ? 

Surely this is a difficult question because its answer would require to 
find the exact solution of the non-relativistic RST eigenvalue problem 
which consists of the (non-relativistic) eigenvalue equations flIV.18a|) - 
flIV.18b|) and the coupled (non-relativistic) Poisson equations flIV.29ap - 



(IIV.29bp . 



B. Hyperfine Splitting 

The effect of level splitting by the magnetic (i.e. spin-spin) interactions 
is experimentally well established and is found to amount to 0, 0008 . . . [eV] 
for the positronium groundstate [8j. It should be clear that such a small 
energy difference between the triplet {^Si) and singlet (^Sq) state falls into 
the order of magnitude of the relativistic effects which such compete with the 
magnetic interaction effects. Therefore it seems very unlikely that the total 
energy difference due to the hyperfine splitting of the groundstate should be 
caused by the magnetic effects alone, but nevertheless it may be interest- 
ing to estimate their relative contribution to the hyperfine splitting of the 
groundstate within the present framework of RST. 

For this purpose, one first has to demonstrate the specific way in which 
this level dichotomy does emerge in RST. Such an effect, however, is rather 
obvious since in the non-relativistic limit one can alternatively put the spin- 
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down component *^^^S^+ 0III.62ap to zero and retain only the spin-up compo- 
nent which itself must then coincide with the first wave amplitude ^^^R+', 
i.e. in place of the former arrangement flIV.58al) - flIV.58bl) one puts now for 
the non- vanishing wave amplitudes 



(V.9a) 



(^)^+ = (2)^+ = . (V.9b) 
This yields the parallelity of both Dirac currents ka{r){= ^""^k^e^j,), 

in contrast to the antiparallelity of the former case (lIV.35aP - (]IV.35bp . Ac- 
cording to the relationships flll.36al) - flll.36bl) and fill. 43aP - fill. 43bp between 



the Dirac currents fca^ and Maxwell currents j"^, the Maxwell equations 
flll.27ap - flll.27bp in three- vector notation flIV.3ip say that the parallelity of 



the Dirac currents ka imply the antiparallelity of the magnetic fields Ha and 
vice versa. These arrangements of the RST fields suggest the following mag- 
netic classification of the positronium states [6j 

h = -h = fcp = ^^^fc^e^ (V.lla) 

ortho-positronium = = k^.J^ = -k2 = k^ (V.llb) 

A,^A, = A^ = ^'^A^e^ (V.llc) 

H,=H^ = H^ = Wif.e; + (')#^e^ , (V.lld) 

and analogously 
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k, = k2 = h = ^'^he^ (V.12a) 

para-positronium = ki = h; j2 = -h = -h (V.12b) 

^'^'^ A = -A2 = ip = (^M^e^ (V.12c) 

= -#2 = Hp = ^"'^HrCr + ^P^H^e^) . (V.12d) 

Thus the present RST ortho-positronium corresponds to the conventional sin- 
glet states (^Sq) and RST para-positronium to the triplet states (^5*1). It is 
true, this RST classification of the positronium states is based upon the (non- 
relativistic) decoupling of the spin-up and spin-down configurations [6] but is 
assumed to hold also for the relativistic case where the spin-up and spin-down 
amphtudes ^"^R^, (")^± remain coupled so that spherically symmetric config- 
urations are not possible, see the eigenvalue equations (lIII.66ap - (]III.66d|) . 



The present magnetic dichotomy of the positronium states lends itself now 
to a very simple calculation of the groundstate hyperfine splitting. Namely, 
for a lowest-order estimate one may resort to the two trial configurations of 
the parallel magnetic-fields flIV.58aP - flIV.58bl) (-^ ortho-positronium) or of 
antiparallel fields (|VJa])-(jVJb]) (-W para-positronium). The corresponding 
magnetic field energy 

E^R^ f UVTb]) differs then only m sign: 



-^iT^lortho = y C^^'^ll-^blP = — -^fT^lpara (V.13a) 

^(™)| _ f ,3^::ft::2_ A(m)| 



-Er Ipara = ~ J ^ ^ll-^pll = ~-^R lortho , (V.lSb) 

provided the magnetic field -f^b/p is computed approximately by means of the 
trial functions for the amplitude combinations {^^^-R+, ^^0+} and {^^^R^, ^'^'>R^}, 
as demonstrated in subsection IV. C {Magnetic Interactions). Thus referring 
to the magnetic interaction energy f ]IV.59p of ortho-positronium, one 
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ends up with the following total energy e!^^: 



Mr'i 



o \ 2 2 



(V.14) 



which of course is the magnetic generalization of the simpler purely electric 
case f lV.Tp . (The upper/lower sign refers to the ortho/para-configurations 
flVmi) -flVl2d|)l 

The minimal value of the total energy E'^\r^) occurs now at the slightly 
shifted position 



Ob 



1+ W1t3 



2«s 



vr 



2aT 



TT 



(V.15) 



and the corresponding minimal value of the energy becomes now in the order 



of cts ^ 



E. 



(0) 

T.min 



4aT 



1 ± 2 (^^ " 

\ TT 



(V.16) 



Therefore the hyperfine splitting IS.e!^'^ is predicted by the present estimate 

as 



'-'-'^T ~ T lortho -'^T I para — R 



n(0)| 



(m) 



2as\' e2 



(V.17) 



TT / 4aB 

This is much smaller than the electrostatic binding energy of 6,80 [eV], cf. 
( IV.Sp . namely 



(0) 



(2, 15 ■ 10"^) ■ 6, 80 [el/] = -1, 46 ■ 10"^ {eV\ . 



(V.18) 



It is true, this is qualitatively in agreement with the experimental fact 
that the binding energy of the ortho-system is greater than that of the para- 
system [8]; but the experimental value of the hyperfine splitting is —8,41 • 
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10~^ [eV] which is six times larger than the present RST prediction (IV.lSp . 
Such a discrepancy may be understood in the sense that for the positronium 
hyperfine sphtting in the order of as ^ it is necessary to use some trial func- 
tion R{r) which is closer to the exact solution than the simple exponential 
function fIV.Qap . In any case, it seems worthwhile to look for the exact so- 
lutions of both the relativistic and non-relativistic RST eigenvalue problem 
in order to test its theoretical accuracy in comparison to the experimental 
situation and the other theoretical approaches. 
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